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FORWARD 


This  report  summarizes  the  results  obtained  to  date  of  an  analysis 
conducted  by  The  University  of  Tennessee  Space  Institute  under  US  Army 
Contract  No.  DAAH01-74-C-0183.  This  contract  was  initiated  under  DA 
Project  No.  IM  26230A  214  and  AMC  Management  Structure  Code  No. 
632303.11.21505.01.  The  technical  effort  was  performed  between  May 
1971  and  September  1975,  under  the  direction  of  the  Aerodynamics  Group, 
Aeroballistics  Directorate,  US  Army  Missile  Command,  Redstone  Arsenal, 
Alabama.  The  Army  technical  representative  was  Dr.  D.  J.  Spring. 
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I. 


INTRODUCTION 


Since  our  last  summary  report  was  prepared  [1],  significant 
progress  has  been  made  both  in  the  areas  of  inviscid  flow  analysis  on 
various  missile  configurations  and  in  understanding  the  viscous  flow  field 
interaction  with  an  exhaust  jet  plume.  Several  publications  have  resulted 
from  the  studies  and  these  are  listed  as  References  2 through  25.  The 
objective  of  the  present  report  is  to  summarize  the  results  of  our  con- 
tinued effort  to  develop  an  understanding  of  the  flow  over  missile  config- 
urations at  transonic  speeds. 

The  goal  of  the  entire  research  project  is  to  develop  calculation 
techniques  that  allow  the  aerodynamic  characteristics  of  missile  configura- 
tions to  be  estimated.  Particular  interest  centers  around  the  effects 
a jet  exhaust  plume  may  have  on  the  flow  field  over  the  body  when  the 
free-stream  conditions  are  transonic.  A strong  interaction  flow  has  to 
be  treated  and  the  basic  flow  configuration  is  shown  in  Figure  1 . 

The  theory  has  developed  within  the  spirit  of  the  component  approach 
wherein  individual  flow  items  are  treated  separately  and  the  results  matched 
together  to  cover  the  entire  flow  field.  The  major  of  these  divisions 
lies  between  the  potential  and  viscous  elements  of  the  flow,  with  each 
being  further  broken  down  to  suit  specific  flow  situations.  Sections  II 
and  III  treat  the  potential  and  viscous  elements  of  the  flow  respectively. 
Section  IV  presents  some  calculated  results  for  plume  induced  separation 
under  various  conditions . 

The  potential  flow  theories  described  in  Section  II  are  based  upon 
the  slender  body  approximation  and  are  restricted  to  small  incidence. 

This  same  restriction  carries  over  to  the  viscous  flow  analysis,  treated 
in  Section  III,  where  the  cross- flow  component  in  the  boundary  layer  is 
neglected. 

Figure  2 is  a summary  statement  of  the  geometric  configurations  that 
are  considered  in  the  potential  flow  analysis.  Each  of  these  items  is 
described  in  greater  detail  in  Section  II  and  their  calculation  in  the 
computer  program  is  indicated  in  Appendix  A.  'The  range  of  applicability 
of  the  theories  associated  with  each  of  these  items  is  expressed  in  Table  1 
as  a function  of  the  pertinent  parameters.  Table  2 presents  similar  infor- 
mation for  the  viscous  flow  calculations.  When  building  up  a configuration 
from  the  various  components,  it  is  the  most  stringent  of  the  parameter 
limitations  which  control  the  region  of  validity  of  the  final  result.  It 
should  be  noted  that  most  of  the  computer  program  components  will  yield 
numbers  when  used  outside  the  recommended  limits  of  Tables  1 and  2 but 
the  results  should  then  be  treated  with  caution. 

Part  II  of  this  report  will  present  the  details  of  the  computer 
program  and  its  utilization. 
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TABLE  1 . PARAMETER  LIMITATIONS  FOR  INVISCID  FLOW 


Component 

M Limits 

a Limits 
(deg) 

Other 

Fore-body 

Ogive 

0.8  - 1.2 

0 to  5 

Very  low  M may  have 

Conical 

0 - 1.5* 

0 to  15 

convergence  problem 
Applicable  till  shock 

Polynomial 

0.7  - 1.5 

0 to  10 

attaches  which  depends 
on  apex  angle* 

Arbitrary  Smooth 

0.7  - 1.5 

0 to  10 

Rear-end 

geometry 

Conical  boattail 

1.0  - 1.5 

0 to  10 

Conical  flare 

0 - 1.5* 

0 to  10 

Same  limitation  as 

Smooth  boattail 

0.8  - 1.5 

0 to  10 

conical  nose* 

Fin 

configurations 
Slender  fin 

0 - 1.5 

0 to  10 

Applicable  to  cruciform 

Nonslender  fin 

0 - 1.0 

0 to  10 

Applicable  to  cruciform 

Control  surface 

Same  as  fin 

C to  30 

Applicable  to  cruciform 

(flap  or  Elevon) 

*Can  be  determined  by  the  graphs  given  in  NACA  - 1134: 
Shock  attachment  angle  and  Mach  number  relationship. 


II.  ANALYSIS  OF  INVISCID  FLOW 

In  order  to  calculate  the  aerodynamic  performance  of  missile 
configurations,  it  is  essential  that  accurate  knowledge  of  the  inviscid 
flow  over  the  body  is  available.  This  is  necessary  as  a starting  point 
in  the  determination  of  the  viscous  flow  over  the  body  and  also  for  the 
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TABLE  2.  LIMITATIONS  FOR  VISCOUS  FLOW  ANALYSIS 


Component 

M Limits 

a,  Limits* 

Other 

Laminar  boundary 
layer 

0 - 2.0 

0 

Unit  Prandtl 
number 

Turbulent  boundary 
layer 

0 ~ 2.0 

0 

Flow  at  slope 
Discontinuity 

- 

0 

Small  angle 
change 

Separation 

>0.5 

0 

Free- shear  layer 

0 ~ 2.0 

0 

Confluence 

> 1.0 

* Approximate  results  can  be  obtained  for  small  incidence 
by  appeal  to  a strip  theory. 


significant  interaction  which  may  take  place  between  the  viscous  and 
potential  flow  components.  The  analysis  starts  with  a discussion  of  the 
inviscid  flow  over  a variety  of  missile  configurations.  The  basic  theore- 
tical work  has  been  properly  documented  in  References  2-9,  11,  12,  15-18, 
and  23  and  will  be  briefly  described  below.  Greater  emphasis  will  be 
placed  upon  a discussion  of  the  flow  over  fin-body  combinations. 

A.  Basic  Potential  Flow  Theories  For  Noses  and  Tails 

The  inviscid,  transonic  flow  small  perturbation  equation  has 
been  employed  to  solve  most  of  the  axisymmetric  nose  and  tail  flows. 
Various  techniques  were  developed  for  solving  the  different  geometries. 

All  theories  developed  so  far  have  been  carefully  compared  with  data  and 
shown  to  possess  very  satisfactory  agreement  with  experimental  data  within 
their  regions  of  validity.  Thus,  all  theories  developed  can  be  employed 
with  confidence  for  performance  predictions  (for  example,  the  Navy  has 
accepted  some  of  our  theories  in  their  standard  program  [26]). 

1.  Flow  Field  Around  Ogive-Cylindrical  Bodies.  The  nonlinear, 
inviscid,  transonic  flow  small  perturbation  equation  has  been  used  to 
solve  the  problem  of  flow  over  an  ogive-cylindrical  body  with  small  angle 
of  incidence.  Hosokawa's  nonlinear  correction  theory  [ 27  and  28]  known 
as  a parabolic  method  with  a stretching  procedure  [2,  4,  8,  16]  has  been 
applied  to  such  geometries.  The  agreement  with  experiment,  in  general,  is  excel 
lent.  This  theory  was  summarized  in  our  last  summary  report  [ 1]  and  a further 
check  with  experiment  was  made  by  Spring  [29] . A typical  comparison  of 
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theory  with  MICOM's  data  is  given  in  Figure  3.  In  these  comparisons, 
the  effect  due  to  the  boundary  layer  build  up  are  also  considered.  From 
this  comparison,  it  is  concluded  that  the  Reynolds  number  effect  on  the 
surface  pressure  calculation  is  negligibly  small  for  the  flow  over  a 
body  of  revolution  at  small  angles  of  incidence  [ 6]  in  the  absence  of  a 
propulsive  jet.  This  is  drastically  different  from  the  situation  that 
can  occur  with  two-dimensional  geometries  [30]. 


2 . Flow  Field  Around  Various  Smooth  Nose-Cylindrical  Bodies. 
The  previously  mentioned  method  has  been  extended  to  cover  general  smooth 
shaped  nose-cylindrical  bodies.  Various  examples  have  been  treated 
[ 16  or  17]  which  range  from  an  arbitrary  smooth  frontbody,  given  in 
tabular  form,  to  bodies  described  by  various  order  polynomials. 


The  inviscid  transonic  flow  small  perturbation  equation  for  a body 
of  revolution  states  that: 


+ \ ^ 
r o 0 


(f>  <t> 

XXX 


(1) 


When  linearized  in  the  parabolic  sense,  the  above  equation  can  be  written 
in  the  form: 


M2)# 
V XX 


+ 


1 _o_ 
r dr 


r de2 


(2) 


The  Hosokawa's  nonlinear  correction  method  assumes  that  the  solution  of 
Equation  1 can  be  obtained  by  considering  the  solution  of  Equation  2 with 
an  addition  of  a correction  function  g,  so  that. 


3>(x,r ,0)  = 0 (x,r,0)  + g(x,r,6) 


(3) 


where  g represents  the  _local  correction.  The  governing  equation  for  the 
correction  term  g,  as  well  as  the  method  of  relaxation  after  passing 
through  the  acceleration  phase  of  the  flow  is  given  in  detail  in  References 
1,  2,  and  8. 


The  necessary  derivatives  of  the  radius  and  the  cross-sectional  area 
of  the  front  body  are  determined  by  a finite  difference  scheme.  The 
m-th  derivative  of  a tabulated  function  is: 


(4) 

I 

Where  D is  the  derivative  operator,  h the  stepsize  between  two  known 
values  of  the  function,  A the  difference  operator,  and  S™  the  Sterling 
number.  The  difference  operator  and  the  Sterling  number  are  respectively; 


Dm  (f (a)  } 


CO 

Z 4 

n = m n 


Anf(a) 


m 


3 
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i (x)  = £(x  + h)  - f(x) 


(5) 


and 


- [£ D"  0 i 


x = 0 


(6) 


The  basic  notation  used  in  the  development  is  shown  schematically 
in  Figure  4.  The  forebodies  studied  are  expressed  by  the  equations: 


R(x)  = D 


(I)2  - d - 1)2] 1/2  - (I 


1/2 


for  a secant  body,  and 
R(x)  = (-I)"  + 1 R 


(I  - r 


+ R 


(7) 


(8) 


for  polynomial  bodies  xvhen  n takes  the  values  3,  4,  or  5.  In  addition  to 
these  configurations,  an  artificially  generated  forebody,  based  upon  a 
tabulated  radius  distribution,  was  used  in  the  computations  [ 17] . 


It  is  found  that  the  sonic  point  (indicated  as  the  parabolic  point 
• in  Figure  4)  as  well  as  the  shock  wave,  shifted  location  considerably  as 
the  forebody  configuration  is  changed  (Figure  5) . The  local  Mach  number 
distributions  for  the  various  bodies  are  shown  in  Figure  6 when  the  free 
stream  Mach  number  is  sonic.  The  effect  of  free  stream  Mach  number  on 
surface  pressure  distribution  for  the  fifth  order  polynomial  body  and  for 
an  ogive  body  are  compared  in  Figure  7. 


3.  Conical  Nose  Cylindrical  Bodies.  An  excellent  inviscid 
flow  solution  has  been  obtained  for  a conically  shaped  forebody  which  can 
be  applied  within  the  subsonic,  sonic  to  supersonic  flow  regions.  The 
originally  proposed  theory  [ 1 and  8]  has  been  modified  [ 15]  to  incorporate 
better  iteration  and  integration  techniques.  The  detailed  analysis  as 
well  as  the  calculative  procedures  are  given  in  Reference  15. 

The  inviscid  transonic  flow  small  perturbation  equation,  as  given 
in  Equation  (1) , is  transformed  by  placing: 

f2(x)  = (l  - M2)  - (7  + 1)M2  0X  . (9) 

Then  it  follows  that: 


f2(x)0  +0  +~0+4r0n=O 

' ' xx  rr  r r 2 66 

r 


(10) 
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where  Equation  (10)  is  a linear  partial  differential  equation  with  a 
variable  coefficient.  This  equation  is  then  solved  to  first  order  appro- 
ximation [15]  with  the  complete  boundary  conditions.  Some  typical  com- 
parisons of  theory  and  experimental  data  are  given  in  Figures  8 and  9. 

Note  that  a blunting  of  the  nose,  as  shown  in  Figure  9,  only  influences 
the  surface  pressure  distribution  near  the  apex.  The  influence  due  to 
the  angle  of  attack  has  also  been  included  in  the  theory  (see  Figures  10 
and  11) . The  theory  allows  the  pressure  distribution  to  be  computed  over 
a rather  wide  Mach  number  range.  The  comparison  with  data  shows  excellent 
agreement  (Figure  12).  A comparison  of  the  changes  in  pressure  drag  with 
Mach  number  is  given  in  Figure  13. 

4.  Flow  Over  a Conical  Flare.  It  has  been  found  that  the 
conical  flare  configuration  can  go  some  way  toward  controlling  the  undesirable 
flow  separation  characteristics  generated  by  a high  thrust  jet.  For 
this  reason  there  is  interest  in  calculating  the  flow  over  this  type  body 
geometry.  The  nose  cone  solution  described  in  the  previous  subsection 
has  been  extended  [ 16  or  18]  . The  agreement  between  the  theory  and 
experimental  data  is  found  to  be  excellent. 

The  treatment  of  the  governing  equation  is  practically  the  same  for 
the  nose-cone  case.  Two  assumptions  are  basic  to  the  development.  First, 
it  is  assumed  that  the  flow  experiences  a stagnation  point  at  the  begin- 
ning of  the  flared  body  and,  second,  that  the  flow  accelerates  to  sonic 
velocity  conditions  at  the  end  of  the  flare.  The  latter  boundary  con- 
dition is  somewhat  verified  by  experiment  [ 31]  , while  the  first  one  is 
only  true  for  inviscid  flow.  In  reality,  the  viscous  flow  will  modify 
the  pressure  field  at  the  flare  root  and  this  is  discussed  later. 

The  conical  flare  geometry  may  be  described  by  the  equation 

R(x)  = (tan  6)  x + RQ  (11) 

where  6 is  the  flare  angle,  Rq  is  the  radius  of  the  cylindrical  body 

ahead  of  the  flare,  and  the  origin  for  the  streamwise  coordinate,  x,  is 
taken  at  the  flare  root.  The  nomenclature  is  shown  in  Figure  14  together 
with  a computed  local  Mach  number  distribution  for  a fixed  flare  angle. 
Comparisons  with  experiment  are  given  in  Figure  15  and  16,  and  the  varia- 
tion of  pressure  coefficient  with  flare  angle  is  given  in  Figure  17. 

5 . Solution  On  a Cylindrical  Body  gehind  a Sharp  Shoulder. 

The  problem  here  is  to  calculate  the  flow  along  a cylindrical  body  just 
downstream  of  the  junction  with  a conical  nose.  Due  to  the  slope  dis- 
continuity at  such  a point  on  the  body,  it  is  necessary  to  employ  dif- 
ferent methods.  The  technique  developed  is  to  combine  the  local  two- 
dimensional  approximation  analysis  proposed  by  Tani  [ 32]  with  the  tran- 
sonic similarity  transformation.  The  detailed  analysis  as  well  as  the 
calculative  procedure  are  given  in  Reference  5,  8,  or  11. 
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The  basic  nonlinear  small  perturbation  equation  for  the  transonic 
flow  over  a body  of  revolution,  Equation  (1),  is  first  transformed  by 
placing 

y = r - R with  y « R 
2 

$ = (1  + ay  + by  + . . .)  9 (x,z)  (12) 

2 ' 3 

z = y + ky  + my  + . . . 


where  cp  is  the  transformed  perturbation  velocity  potential  and  a,  b,  c, 

. . . , k,  m,  n,  . . . are  constants.  These  constants  can  be  chosen  so 
that  Equation  (1)  can  be  transformed  into  the  corresponding  two-dimensional 
form 

(l  - M2)  cp  + cp  = (1  + y)M2  cp  cp  . (13) 

The  body  boundary  condition  takes  the  form 

- § 9 (x,0)  + 9Z  (x,0)  = H . (14) 

Equation  (13)  is  now  solved  for  the  case  of  a sonic  flow,  and  a wider 
, transonic  flow  range  obtained  by  applying  Spreiter's  similarity  transfor- 
mation [33],  A typical  comparison  of  this  theory  with  experimental  data 
is  given  in  Figure  18.  The  agreement  is  very  good  in  the  upper  transonic 
flow  regime.  In  the  lower  transonic  flow  regime  it  is  found  that  a shock 
wave  and  a local  flow  separation  are  usually  present.  The  location  of 
the  shock  wave  may  be  estimated,  therefore,  by  employing  a turbulent 
boundary  layer  separation  criterion  to  recover  the  free  stream  pressure 
[5  and  34] . The  result  of  employing  such  a scheme  is  compared  with  the 
experiment  on  Figure  19. 

6.  Flow  Over  a Conical  Boattail.  The  above  analysis  is  also 
applicable  to  a conical  boattail  geometry  as  long  as  the  approaching  Mach 
number  is  close  to  unity.  To  develop  the  solution  for  a conical  bottail, 
the  body  geometry  is  expressed  in  the  form 

R(x)  = 1 - Ax  (15) 

where  A denotes  the  boattail  angle  in  radians. 

Now  the  analysis  for  a body  with  discontinuous  curvature  can  be 
applied  directly.  It  is  necessary,  however,  to  assume  that  the  body  is 
sufficiently  long  for  the  flow  ahead  of  the  boattail  junction  to  be 
considered  as  a uniform  free  stream.  The  detailed  analysis  is  referred 
to  in  References  5 and  8.  The  calculation  terminates  when  the  free- 
stream  Mach  number  is  recovered  if  this  condition  is  satisfied  before 
the  end  of  the  body  is  reached. 
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A typical  comparison  with  MICOM's  data  is  shown  in  Figure  20  where 
it  is  seen  that  the  agreement  is  good. 

7.  Flow  Over  a Smoothly  Contoured  Boattail.  The  above  analysis 
has  been  extended  to  a smooth  boattail  by  simply  assuming  that  the  boat- 
tail  may  be  segmented  into  a series  of  conical  boattails  each  at  the  local 
slope.  This  analysis  was  carried  out  by  Whitfield  [35]  and  a typical 
calculation  is  shown  in  Figure  21. 

B.  Fin-Body  Interaction  By  Slender  Body  Theory 

1,  Slender  Body  Theory  and  Flow  Model.  The  research  into 
fin-body  interaction  problems  has  assumed  increased  importance.  In 
response  to  this  need,  an  analysis  of  the  aerodynamic  performance  of  a 
slender  fin-body  combination  has  been  developed.  This  permits  us  to 
estimate  the  pressure  coefficient  on  slender,  delta-cruciform  fins  and  on 
cylindrical  bodies  in  inviscid  flow.  The  analysis  has  been  extended  to 
include  the  case  with  small  angle  of  attack  and  with  sideslip  (due  to 
yaw)  as  well  as  the  case  of  a canted  fin.  The  analysis  is  mainly  based 
on  the  slender  body  theory  of  Spreiter  [ 36  and  37]  and  the  source  and 
sink  method  of  Adams  and  Dugan  [38] . 

In  the  light  of  slender-body  theory,  the  current  analysis  (for  the 
details,  see [23])  is  valid  not  only  for  subsonic  flow  but  also  for  mod- 
erately supersonic  flow  over  a fin-body  combination;  provided  that  the 
fin  has  a subsonic  leading  edge.  Furthermore,  the  analysis  is  valid 
even  at  free-stream  speeds,  very  close  to  sonic,  provided  that  no  shock 
waves  appear  on  the  fin-body  combination.  The  analysis  developed  here 
checks  well  with  experimental  data  over  the  Mach  number  range  mentioned 
and  details  will  be  given  later.  The  analysis  is  outlined  below. 

Much  previous  research  on  fin-body  combination  problems  has  centered 
around  finding  the  "load"  distribution  and  little  work  is  available  con- 
cerning the  surface  pressure  distribution.  Spreiter  [36  and  37]  developed 
a method  based  upon  slender  body  theory  for  calculating  the  aerodynamic 
loading  on  slender  delta  fin-body  combinations  with  small  angles  of  attack 
and  sideslip  (due  to  yaw).  Adams  and  Dugan  [38]  have  treated  the  load 
distribution  on  a slender  canted  fin  by  means  of  the  source  and  sink 
method.  Nielsen  [39]  summarized  previous  work  and  considered  the  aero- 
dynamic loading  in  a more  general  case.  He  treated  slender  fin-body 
combinations  with  small  angles  of  attack  and  sideslip  (due  to  yaw)  in 
addition  to  fins  with  small  cant  angles.  Recently,  Spreiter  and  Stahara 
[40,  41]  investigated  in  detail  the  pressure ’coefficient  distribution  on 
slender  fin-body  combinations  at  transonic  speed.  However,  for  the  case 
of  the  flow  over  a slender,  zero  thickness  fin  mounted  on  a cylindrical 
body  reduces  to  Spreiter 's  earlier  works  [36,  37]  and  to  that  of  Lawrence 
and  Flax  [42] . 
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Figure  22  shows  three  different  body- fixed  Cartesian  coordinate 
systems  used  in  the  present  analysis.  The  Cartesian  coordinate  system 
with  the  axes  x,  y,  and  z is  taken  such  that  the  x-axis  is  positive  rear- 
ward and  coincides  with  the  longitudinal  fin-body  combination  axis,  the 
y-axis  is  positive  to  the  right,  facing  forward,  and  can  rotate  about  the 
x-axis  lying  in  a pair  of  horizontal  fins  of  symmetry,  and  the  z-axis  is 
positive  vertically  upward  and  can  rotate  about  the  x-axis  lying  in  a 
pair  of  vertical  fins  of  symmetry.  The  x,  y,  ~z  coordinate  system  is 
taken  where  the  x-axis  aligns  parallel  to  the  free  stream  velocity,  the 
"z-axis  is  in  the  plane  containing  the  x-axis  and  the  fin-body  longitudinal 
axis.  The  Cartesian  x1 , y',  z1  coordinate  system  is  taken  where  the  x'- 
axis  is  positive  rearward  of  the  body  and  coincides  with  the  fin-body 
longitudinal  axis,  and  the  x'-axis  is  positive  upward  normal  to  the  body. 

The  angle  otQ  is  the  angle  between  the  x'-axis  and  the  x-axis.  The 

angle  cp  is  the  angle  of  bank  and  can  rotate  about  the  x'-axis.  For  small 
ac,  the  angle  of  attack,  ot,  and  sideslip  (due  to  yaw)  p,  can  be  defined 

as  follows: 

a = ac  cos  cp  ) 

> ' (16) 

P = ac  sin  cp  \ 

The  free  stream  velocity,  VQ , can  be  decomposed  into  two  components 

for  small  q;^;  the  axial  flow  with  the  velocity  VQ  and  the  cross  flow, 

V a . Furthermore,  the  cross  flow,  V a , can  be  decomposed  into  two  com- 
o c o c 

ponents  with  velocities  VQa  and  VQP,  respectively.  In  the  case  of  a fin- 

body  combination  where  the  fin  has  a small  cant  angle  about  the  hinge  axis 
at  the  mid  point  of  the  root  chord  length,  the  cant  angle  will  be  desig- 
nated by  5 as  shown  in  Figure  23. 

A positive  rolling  moment  contributed  to  a left  fin  with  its  trailing 
edge  inclinded  downwards  will  be  designated  as  positive  "cant"  in  this 
analysis . 

Figures  24  and  25  show  the  computational  flow  models  and  coordinate 
systems  used  in  Sections  2 and  3 respectively.  In  these  figures,  "a"  is 
the  body  radius  and  "s"  the  fin  semispanwise  length.  Figure  26  shows  the 
computational  model  with  a fin  having  different  ratios  of  semispan  to  body 
radius  used  in  Section  4. 

The  perturbation  velocity  potential  equation,  which  governs  the  flow 
around  the  slender  fin-body  combination,  can  be  written  as 

0 +0  =0  (17) 

yy  zz  ' ' 
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in  the  x,  y,  z coordinate  system  [39] , where  0 is  the  small  perturbation 
velocity  potential,  and  the  subscripts  indicate  partial  derivatives  with 
respect  to  the  corresponding  coordinate  axes.  This  means  that  the  equa- 
tion which  governs  a flow  around  a slender  body  remains  invariant  for  a 
slight  rotation  of  the  coordinate  system.  The  boundary  conditions  for 
0 in  the  fin-body  combination  can  be  written  as 


0 (°°)  = 0 , 
in  the  far  field,  and 


(18a) 


n 


2 


y'j  + (V„« 


+ 0 


Z) 


= 0 


(18b) 


on  the  body.  Here  n^  and  n^  are  the  directional  cosines  of  the  unit 
vector  normal  to  the  body  surface  in  the  y and  z directions,  respectively 


0 = -V  a ± V 5 

zoo 


(18c) 


on  the  horizontal  fins,  where  the  upper  sign  is  for  the  right  horizontal 
fin  and  the  lower  sign  is  for  the  left  horizontal  fin,  and 

0 = V p T V 6 (18d) 

y oK  o 

on  the  vertical  fin  where  the  upper  sign  is  for  the  upper  vertical  fin 
and  lower  sign  is  for  the  lower  vertical  fin,  respectively. 

The  perturbation  potential  for  a unit  velocity  can  be  decomposed  as 


0 • a0a  + e0p  - 60s 


(19) 


where  0 , 0_ , and  0 are  the  perturbation  potentials  for  unit  angle.  The 

Ctp  o 

expression  for  0 and  0 can  be  written  as 
ct  p 


a 


- 2 + 4a4  cos2  2 0 


+ 


+ 


2 

s 


4 

s 

% 


(20) 
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where  the  upper  sign  is  for  the  upper  half  of  a fin-body  combination  as 
shown  in  Figure  27  (i.e.,  0 < 0 < tc)  , and  the  lower  sign  is  for  the  lower 
half  (i.e.,  jt  < 0 < 2ft).  For  0 the  expression  is 

/ - P 


0 


0 


2 

r 


cos  2 0 


2 

cos 


20 


+ s 


+ 


M) 

P4) 

2 

r cos  2 0 

1/2 

V B / 

V r / 

- 

1/2 


+ y 


(21) 


where  the  upper  sign  is  for  the  right  half  of  a fin-body  combination 

TC  tf 

(i.e.,  - ^ < 0 < , and  the  lower  sign  is  for  the  left  half  (i.e., 

f ^ 0 < | Tt)  . 


By  a conformal  transformation,  the  flow  in  the  physical  plane  (termed 
the  x-plane)  is  mapped  into  a flow  in  the  mapping  plane  (termed  the  a- 
plane)  (Figure  28)  by  using  a source  and  sink  method  [38]  , the  perturba- 
* tion  velocity  potential  0 for  unit  velocity  and  unit  cant  angle,  can  be 
written  as 


*8 


sin  2 (7  - p) 
sin  2 (7  + p) 


+ *2  + ^3  + ^4) 


sin  2 p 


(22) 


for  all  values  of  p,  (-jt  < p < jt)  . In  this  equation  and  are  the 

complete  elliptical  integrals  of  the  third  kind,  while  and  1^  are  the 

incomplete  elliptical  integrals  of  the  third  kind.  (For  details  of  the 
theory,  see  [23]). 


The  pressure  coefficient  C^,  including  quadratic  terms  but  neglecting 

the  higher  order  terms  in, the  perturbation  velocity  components,  can  be 
written  for  a unit  free  stream  velocity  in  the  form: 


C = -2  (u 
P 


2 2 

- vf3  + wa)  - (v  + w ) 


(23) 


where  u,  v,  and  w are  the  perturbation  velocity  components  parallel  to 
the  x,  y,  and  z axes,  respectively. 
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The  load  distribution  on  the  fin-body  combination  in  this  analysis 
is  defined  as  follows: 


Ap/q  (^p)lower  (^p)upper 
(Cp)left  " (^p) right 


(24) 


2‘, 


where  q is-  the  dynamic  pressure  (equal  to  PqVq  / 2),  and  the  subscripts 

"lower"  and  "upper"  are  t)ie  lower  (negative  z)  and  upper  (positive  z) 
part  of  the  fin-body  combination.  The  subscripts  "left"  and  "right" 
are  the  left  (negative  y)  and  right  (positive  y)  part  of  the  fin-body 
combination,  respectively,  when  facing  forward  (see  Figure  24a).  By 
considering  the  symmetric  properties  of  the  perturbation  velocity  compo- 
nents on  the  upper  and  lower  surfaces,  it  is  found  that  the  load  distribu- 
tion on  the  horizontal  fin  is  given  as 


-f  ' + 4* 

Ap/q  = -4a  u^  + 4S  Ug  + 4a  pv^ 


(i  - V)  i1 ■ V)  v (25) 


where  the  "+"  superscript  on  the  perturbation  velocity  components  refers 
to  the  lower  panels  of  the  horizontal  fins. 

* I 

The  load  distribution  on  the  vertical  fins  is  given  as 

+> 


+ + + 
Ap/q  = -4p  u + 46  u.  - 4a  wQ 

pop 


(l  + wa+)  + ta&  w5+  (l  + wa+)  (26) 


where  the  "+"  superscript  on  the  perturbation  velocity  components  refer 
to  the  right-hand  panels  of  the  vertical  fins. 


The  load  distribution  on  the  body  (positive  when  upward)  is  expressed 
in  the  form 

Ap/q  - -ta  ua+  + 46  UB+  - tap  jy-  (l  + y)  - va+  (l  - vp+ 

- 4P6  [^v6+  (l  - vp+)  - W6+  y ], 

where  now  the  "+"  superscript  on  the  perturbation  velocity  components 
refer  to  the  lower  part  of  the  body. 

2 . Surface  Pressure  Distributions  On  Fin  Body  Combinations 
With  Angles  of  Attack  and  Yaw. 


(a)  Effects  of  Angles  of  Attack  (with  an  example  for 
a = 3°,  p - -5  = 0°). 
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An  example  o£  the  computation  of  the  surface  pressure  and  the  load 
distribution  was  carried  out  on  a fin-body  combination  at  angle  of  attack 
with  zero  fin  cant  angle.  The  flow  model  used  is  shown  in  Figure  24a. 
Figure  29  presents  the  calculated  C distribution  corresponding  to  each 
side  of  the  fin  panel.  P 

The  C distribution  pattern  is  symmetric  which  respect  to  the  z-axis 

which  aligns  in  the  cross  flow  direction.  The  expected  reason  why  the 
positive  C (i.e.,  the  pressure  side)  close  to  the  horizontal  fin  edge 
P 2 

decreases  is  that  the  second  order  term  v in  Equation  (23)  is  getting 
larger  when  compared  to  the  first  order  velocity  components  near  the  fin 

edge.  The  C distribution  on  the  vertical  fins  is  the  same  and  will 
P 

become  zero  as  z becomes  ± °°.  Figure  30  shows  the  body  distribution 

which  is  also  symmetric  with  respect  to  the  z-axis.  It  shows  the  greatest 
change  at  the  junction  of  the  horizontal  fins  with  the  body.  The  schematic 
load  distribution  on  the  entire  fin-body  combination  is  shown  in  Figure  31. 
The  body  has  an  ogive  nose  in  this  example.  It  should  be  noted  that  no 
attempt  was  made  to  calculate  the  load  distribution  for  the  vertical  fins 
in  this  particular  case.  Thus,  the  load  distribution  on  the  cruciform 
fin  is  exactly  the  same  as  that  for  the  planar  fin  case. 

(b)  The  Effects  of  Angles  of  Yaw  (with  an  example  of 
p = 3°,  a = -5  = 0°): 

Again,  as  a demonstration  of  the  calculation  procedure,  the  sideslip 
angle  was  taken  to  be  p = 3°,  while  OL  = -S  = 0°  . In  this  case,  the 

distribution  and  the  load  distribution  will  have  exactly  the  same  pattern 
as  would  be  obtained  by  rotating  the  C and  the  load  distribution  patterns 

of  the  above  example  (see  Figures  29,  30,  and  31)  90°  counterclockwise. 

(c)  The  Combined  Effects  of  Angles  of  Attack  and  Sideslip 
(with  an  example,  a = p = 3°,  -5  = 0°): 

As  a sample  calculation,  the  values  a.  = P = 3°  were  assumed.  In 
this  simple  combined  case,  it  is  expected  that  the  distribution  and 

load  distribution  will  have  approximately  the  same  pattern  as  the  resultant 
pattern  of  the  above  cases  for  VQa  and  VQp. 


Figures  32  and  33  show  the  C distributions  on  the  surfaces  of  the 

P 

fins  and  the  body,  respectively.  The  flow  model  for  this  case  is  shown 
in  Figure  24b. 


The  0^  distribution  is  symmetric  with  respect  to  the  z'-axis  which 

has  the  same  direction  as  the  resultant  cross  flow  V a with  bank  angles 

o c 

of  45°,  because  a = p.  The  reason  that  the  C distribution  on  the  surface 
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has  a larger  negative  value  (i.e.,  on  the  suction  side)  and  larger  posi- 
tive <1!^  (i.e.,  on  the  pressure  side)  in  the  second  quadrant  and  the  fourth 

quadrant,  respectively,  as  shown  in  Figure  32  and  33,  is  because  the 
pressure  caused  by  the  cross  flows  V^a  and  intensify  each  other  in 
these  quadrants. 

By  examining  the  distribution  on  the  body,  an  interesting  charac- 
teristic is  found  (Figure  33) . The  value  of  at  one  point  on  the  first 

or  the  third  quadrant  of  the  body  becomes  "zero".  This  means  that  the 
pressure  at  these  points  recovers  to  that  of  the  free  stream  and  is  a 
consequence  of  the  co-existence  of  both  compression  and  suction  regions 
in  the  same  quadrant.  Figure  33  presents  the  quantitative  surface  pressure 
distributions  around  the  circumference  of  the  body. 


3.  Effect  of  Canted  Fin. 

(a)  C and  Load  Distribution  on  Cruciform  Canted  Fin-Body 

P 

Combination  (with  an  example  for  -5  = 3°,  a = (3  = 0°): 

As  a sample  calculation,  the  pressure  and  load  distribution  were 
calculated  for  cruciform  fins  at  a cant  angle  of  "3  " (in  a negative  sense, 
see  Figure  25)  on  a circular  cylinder  with  its  longitudinal  axis  aligned 
to  the  free  stream  direction.  The  geometry  is  shown  in  Figure  25a.  The 
computed  C^  distribution  on  the  fins  is  shorn  in  Figure  34.  The  C^  dis- 
tribution on  the  corresponding  panel  of  each  fin  has  exactly  the  same 
pattern  and  has  a symmetric  pattern  with  respect  to  the  origin  of  the 
coordinate. 

The  Cp  distribution  on  each  quadrant  of  the  body  also  has  the  same 

pattern  which  is  again  symmetric  with  respect  to  the  origin  of  the  coor- 
dinates. The  Cp  distribution  appearing  in  Figure  35  is  regarded  as  the 

perturbation  disturbance  created  by  the  fins  since  the  body  alone  would 
not  generate  a disturbance  The  reason  why  the  Cp  distribution  on  each 

quadrant  has  a zero  point  is  a consequence  of  the  pressure  differential 
across  the  fins  (see  Figure  34) . . 

(b)  Comparison  with  the  Case  When  Only  the  Horizontal 
Fins  are  Canted  (with  an  example  for  -6  = 3° , a = 0 = 0°): 

Figure  36  shows  the  computed  C distribution  on  the  fins.  The  C 

P P 

distribution  shows  nearly  the  same  values  as  those  in  Figure  29  of 
subsection  (a) . 
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Figure  37  shows  the  distribution  on  the  body  which  results  from 
the  pressure  disturbance  transmitted  by  the  horizontal  fins.  The  dis- 
tribution on  each  quadrant  (Figures  30  and  37)  yields  a similar  pattern. 
The  approximate  distribution  on  both  the  fins  and  the  body,  in  the 

case  where  all  four  fins  are  canted,  can  be  obtained  by  superimposing  the 

C distribution  for  the  cases  w 

P 

are  canted  (Figures  36  and  37)  . 


0^  distribution  for  the  cases  when  only  the  vertical  and  horizontal  fins 


(c)  Comparison  of  C and  Load  Distribution  on  Fin  Body 


Combinations  in  Various  Cases: 


The  Cp  and  load  distributions  for  various  cases  have  been  compared. 

For  one  of  these  examples,  Figure  38  shows  the  C^  distribution  on  the 

right  horizontal  fin  and  body  while  Figure  39  shows  the  load  distribu- 
tion on  the  fin-body  combinations.  The  load  distribution  on  the  body  is 
computed  for  a unit  projected  spanwise  length.  The  reason  why  Ap/q  on 
the  body,  for  the  case  of  all  fins  canted,  is  equal  to  zero  at  y/a»  0.7 
is  the  distribution  in  this  case  becomes  the  same  for  both  the  upper 

and  lower  parts  of  the  body  (see  Figures  35  and  38)  at  the  midpoint  of  a 
half  body  (y/a  = cos  45°  « 0.71). 

4.  Comparisons  Of  Various  Configurations  Of  Fin-Body 
Combinations , 


2(a)  Surface  Pressure  Distributions: 


As  a sample  calculation,  the  case  when  a = g = -8  = 3 has  been  con- 
sidered. It  can  be  surmised  that  the  C and  Ap/q  distribution  in  this 

P 

case  will  yield  approximately  the  same  pattern  as  those  obtained  from  the 
results  of  the  combined  cases  of  a = (3  = 3°  (Section  c ) and  -5  = 3° 
(Section  a ).  In  fact,  the  computed  C^  distribution  as  shown  in  Figures 

40  through  42  is  nearly  the  same  as  that  of  the  combined  C distribution 

for  the  case  a = p = 3°  (Figures  32  or  33),  and  for  the  case  -5  = 3° 
(Figures  34  or  35) . The  C^  distributions  on  each  surface  are  more  similar 

than  those  of  the  corresponding  case  for  a,  = (3  = 3 due  to  the  small 
influence  of  the  -6  = 3°  cant  angle  in  such  cases.  A larger  value  of 
positive  0^  is  found  in  the  fourth  quadrant  while  the  negative  C^  is 

larger  in  the  second. 


2(c)  Mutual  Effects  of  a,  (3,  and  8 — Discussion  on 

Coupling  Terms: 

There  will  be  a "mutual  interaction"  between  the  effects  of  angle  of 
attack  oc,  the  angle  of  yaw  |3,  and  the  cant  angle  5.  The  magnitude  of  this 
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interaction  flow  field  has  been  estimated.  Figures  43  and  44  show  the 
load  distribution  caused  by  each  perturbation  angle  for  one  particular 
configuration  at  s/a  = 2.5.  These  figures  also  show  the  interaction 
between  the  three  effects  of  incidence , yaw  and  cant  angle.  From  these 
figures,  it  is  observed  that  the  major  contribution  to  Ap/q  on  each  fin 
is  made  by  Ap/q  on  the  right  and  left  horizontal  fins  from  the  incidence 
effect,  and  on  the  upper  and  lower  vertical  fins  from  the  yaw  effect. 

This  is  mainly  due  to  the  first  term  on  the  right-hand  side  of  Equations 
(25)  and  (26),  respectively.  In  general,  the  effect  of  cant  angle  on 
Ap/q  is  small,  and  is  due  to  the  second  term  on  the  right-hand  side  of 
Equations  (25)  or  (26) . Pertaining  to  the  load  distribution,  the  inter- 
action between  incidence  and  yaw  affects  all  fins,  but  the  mutual  inter- 
action between  incidence  and  cant  angle  appears  only  on  the  vertical  fins. 
The  mutual  interaction  between  yaw  and  cant  angle  appears  on  the  hori- 
zontal fins  only,  as  can  be  realized  by  examining  Equations  (25)  or  (26) . 
It  is  clear  that  the  coupling  terms  mentioned  above  do  not  contribute 
significantly  to  the  load  distribution  over  the  fins. 

C.  Fin-Body  Interaction  By  Singularity  Distribution 

1.  Basic  Considerations.  The  aerodynamic  performance  of 
various  kinds  of  fin-body  combinations  in  inviscid  flow  can  be  analyzed 
by  the  "singularity  method"  [43  and  44] . Korner's  approach  [43]  is 
followed  in  this  analysis.  In  general  the  singularity  method  has  advan- 
tages over  the  slender  body  theory  since  there  is  less  restriction  on 
admissible  fin  shapes.  The  pressure  coefficient  can  also  be  calculated 
rather  simply. 

The  surface  pressure  distribution  on  rectangular,  sweptback,  and 
tapered  sweptback  fin-body  combinations  have  been  obtained.  The  analysis 
has  also  been  extended  to  include  slender,  cruciform  and  canted  delta 
fin-body  combinations  as  discussed  earlier,  and  also  to  include  the 
case  of  slender,  cruciform,  clipped  delta  fin-body  combinations  with  a 
deflected  small  control  fin.  The  results  have  been  compared  with  slender 
body  theory  and  with  other  singularity  methods  [44] . 

The  present  method,  with  slight  modification,  can  also  be  applied 
to  a fin-body  combination  when  the  fin  is  twisted,  cambered,  of  finite 
thickness  and  of  arbitrary  plan  form. 

For  simplicity,  the  following  assumptions  were  made  in  the  present 
analysis,  i.e.,: 

(i)  The  free- stream  flow  is  subsonic,  steady,  inviscid  and 
uni form. 

(ii)  The  fin  is  of  zero  thickness  and  has  straight  leading  and  , 
trailing  edges  and  a streamwise  tip. 
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(iii)  The  fins  are  attached  to  a cylindrical  body. 

(iv)  The  angle  of  attack  is  small. 

Based  on  the  singularity  method  concepts,  constant  strength  vortex 
singularities  on  a horizontal  plane  are  used  to  replace  the  lifting  fin. 
Image  vortices  are  placed  inside  the  body  in  order  to  compensate  the  body 
displacement  effect.  To  obtain  a good  aerodynamic  interaction  for  a fin- 
body  geometry,  source  — sink  singularities  are  distributed  along  the  body 
surface  and  doublet  singularities  are  distributed  along  the  body  axis. 

The  doublet  distribution  gives  the  effect  of  the  body  being  inclined  to 
the  free  stream.  An  interaction  scheme  is  used  so  that  all  the  simulated 
lifting  surfaces  can  eventually  satisfy  the  flow  boundary  conditions 
everywhere  on  the  surface  of  the  fin-body  combination. 

Figure  45  shows  a schematic  model  of  a basic  fin-body  combination 
for  this  analysis,  A fin  (or  wing)  has  been  divided  into  many  panels 
and  each  one  has  been  replaced  by  a horseshoe- type  vortex  which  consists 
of  three  vortices  (i.e.,  one  bound  vortex  and  two  free  vortices).  The 
body  has  also  been  divided  into  many  panels,  each  consisting  of  a correct 
part  of  the  body,  and  it  has  been  replaced  by  a source  or  sink  on  its 
surface. 

2.  Basic  Equations  and  Coordinates.  Figure  46  shows  the 
coordinates  used  on  the  fin.  Each  fin  is  divided  into  a number  of  trape- 
zoidal panels  having  their  sides  parallel  to  the  free  stream.  Each  panel 
has  three  vortices  of  constant  strength,  i.e.,  one  bound  vortex,  which 
is  located  at  the  one  quarter  local  panel  chord  length,  and  two  free 
trailing  vortices  which  are  located  at  the  local  panel  edges  and  are 
parallel  to  the  free  stream.  The  control  point  location,  at  which  the 
boundary  condition  is  to  be  satisfied,  is  assumed  to  be  the  three  quarter 
point  of  the  local  panel  chord  length  and  on  its  center  line,  shown  in 
Figure  47.  In  addition  to  this,  the  Goethert  compressible  similarity 
rule  has  been  employed  in  the  present  analysis,  i.e.,  the  coordinates 
normal  to  the  free  stream  have  been  multiplied  by  the  Prandtl-Glauert 
parameter  (3.  The  following  five  steps  are  essential  to  the  analysis  and 
were  established  in  [43]: 

(i)  Determining  the  circulation  distribution  (7^)  on  the  fin: 

The  normalized  perturbation  velocity  components  at  any  control 
point  on  panel  "i"  induced  by  any  panel  "j"  (that  is  four  discrete 
horseshoe  vortices,  two  from  a pair  of  fins  themselves,  and  two  from 
their  body  images)  can  be  written  as, 


19 


u. 

n 

1 

U 

CO 

E 

J=i 

V, 

n 

i 

U 

CO 

E 

j=i 

W. 

1 

U 

00 

n 

E 

j=l 

3-J  j 


‘ij  'j 


_1  IJ  J 


i = 1,  n 


where:  u^,  v^,  w^  are  the  perturbation  velocity  components  and  P^, 

Q..  and  R. . are-  the  aerodynamic  matrices, 
ij  iJ 

I 

The  third  equation  as  stated  above  is  used  for  satisfying  the  boundary 
conditions  on  any  panel,  i.e., 


w.  n 

-Bio.  = ~ = E R.  .7 . 
i U ij  1 

CO  J = 1 J J 


i = 1,  n 


where  = a±±’  Equation  (29)  consists  of  n simultaneous  linear  equations 

with  n x n constant  coefficients.  The  7.  can  be  solved  by  a Gaussian 
elimination  technique . ^ 

(ii)  Determining  the  source  distribution  (qy)  on  the  body 

The  induced  velocity  components  at  any  point  on  panel  "v"  on  the 
body  surface  due  to  the  fin  and  its  image  vortices  associated  with  the 
control  point  of  panel  "j"  is  expressible  as. 


% ’ vbv  - £ Evi  - - 1.  2m 


In  cylindrical  coordinates  there  exist  the  transformation: 

v . = v,  cos  0 + w,  sin  0 
nl  b b 


v.  = -v,  sin  0 + w,  cos  0 

t D D 


where  the  geometrical  relationships  are  illustrated  in  Figures  48  and  49. 
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The  m-th  approximation  of  the  source  strength  q^  at  a point  v on 

the  body  surface,  which  cancels  the  induced  normal  velocity  v ^ in 
Equation  (31) , can  be  written  as  l' 


211  q / u„ 

V " -2(Vnl)r»H  - y£  X 
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and 


as 


A Sy  = p a A0  Axy 


(33) 


(iii)  Determining  the  downwash  (w,  ) on  a fin  due  to  sources 

D p 


distributed  on  the  body  surf dee 

The  downwash  on  panel  "p"  due  to  source  on  panel  "v"  can  be  written 
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(34) 


(iv)  Determining  the  circulation  distribution  (7.)  for  the  second 
iteration 


As  a result  of  the  downwash  obtained  in  Equation  (34),  the  corre- 
sponding angle  of  attack  induced  on  a fin  can  be  written  as  w^/U^,  This 

term  should  be  added  to  the  geometrical  angle  of  attack  in  Equation  (29) 

so  that  co.  is  written  as 
x 


“Wi  " -°fi  + U J 1 " n (35) 

00 

After  these  co^,  values  have  been  obtained  on  the  fin  panels,  the 
procedure  is  repeated  until  the  values  of  7^  do  not  change  greatly. 

(v)  Extension  to  the  case  of  a body  with  angle  of  attack. 
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Downwas(h  on  a fin  due  to  the  doublet  distribution  along  the  body 
axis  (as  shown  in  Figure  50)  can  be  written  as 
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(36) 


This  downwash  has  been  assumed  constant  at  any  spanwise  location,  where 
it  causes  an  apparent  change  in  the  angle  of  attack  (un)  as  mentioned 

above.  Thus,  in  this  case,  ah  should  be  replaced  by, 
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The  aerodynamic  performance  coefficient  can  be  computed  in  the  usual 
manner  for  the  pressure  coefficient,  C on  the  panel, 


for  fin 


for  body 


(38) 


where  "+"  denotes  the  upper  surface  of  the  fin  and  the  lower  surface. 


For  load  distribution,  AC  , there  results 
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The  spanwise  lift  distribution,  CT (y)  is  given  as 


cL(y) 


2b 

P2i(y) 


(40) 


where  n'  is  the  total  number  of  panels  along  the  chord  l(y)  at  a fixed 
y station. 
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For  the  lift  coefficient,  C.^, 
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the  final  result  is 
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and  represents  the  projected  area  on  the  horizontal  plane. 

3 . Results  and  Discussion  on  Various  Fin-Body  Combinations. 


(a)  Rectangular  Fin-Body  Combination: 

The  spanwise  lift  distributions  have  been  calculated  for  a rectangu- 
lar fin  alone  and  for  a fin-body  combination  configuration  at  =0. 

They  agree  well  with  the  original  calculations  of  Korner  [43]  as  is  shown 
in  Figure  51.  The  discussion  of  the  compressibility  effect  will  be  given 
later.  The  body  treated  in  this  example  has  zero  angle  of  incidence,  and 
the  fins  are  canted  by  1 (clearly,  results  for  a larger  cant  angle  can 
be  obtained  by  a simple  linear  multiplication  of  the  amplitude).  The 
computed  results  for  the  fin-body  interaction  agreed  well  with  the 
limited  experimental  data  shown  in  Figure  51,  except  in  the  region  very 
close  to  the  fin  root.  Korner  pointed  out  that  this  disagreement  resulted 
from  the  presence  of  the  boundary  layer  in  the  real  flow  over  the  fin  and 
body  [ 43] . This  suggests  that  a study  of  viscous  effects  is  necessary  in 
order  to  obtain  a more  accurate  result  in  this-  region.  As  can  be  seen 
from  Figure  51,  when  comparing  cases  with  and  without  body,  the  mutual 
interaction  between  the  body  and  the  fin  is  rather  strpng. 


The  chordwise  load  distributions  on  the  body,  for  the  same  fin-body 
geometry  as  shown  in  Figure  51,  were  computed  at  three  different  azimuthal 
positions.  The  results  are  shown  in  Figure  52.  It  can  be  seen  that  the 
influence  of  the  fin  on  the  body  becomes  stronger  near  the  fin  root.  It 
affects  both  the  upstream  and  downstream  flow  beyond  the  fin  location. 

The  calculated  chordwise  lift  distribution  on  the  fin  very  close  to  the 
fin  root  is  also  included  in  Figure  52  for  a comparison. 

(b)  Constant  Chord  Sweptback  Fin-Body  Combination: 

Figure  53  shows  a calculated  chordwise  load  distribution  on  the  fin 
of  a 45  sweptback  fin-body  combination.  The  variation  of  the  load  dis- 
tribution, which  is  rather  small  along  the  spanwise  direction,  can  be 
seen  from  this  figure. 
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(c)  Tapered,  Sweptback  Fin-Body  Combination: 

The  sweptback  fin-body  combination,  with  two  different  taper  ratio 
fins,  as  used  by  Korner  [43],  has  been  investigated.  Figure  54  shows  the 
computed  spanwise  lift  distribution  on  this  type  of  geometry  when  the 
taper  ratio  X = la/li  is  0 and  one-third  (where  1 li ' is  the  maximum 
chord  length  through  the  body,  and  'la'  is  the  fin  tip  chord  length). 

The  sweptback  angle  for  both  cases  was  taken  as  30°  at  one-quarter  of 
the  local  chord  length.  The  angle  of  incidence  of  the  body  was  assumed 
to  be  1°  (again,  for  larger  angles,  a linear  multiplication  of  the  amp- 
litude is  sufficient)  and  the  fin  was  assumed  to  have  a zero  cant  angle. 
The  lift  distribution  for  X = 0 agrees  well  with  the  characteristics  of 
the  delta  fin-body  combination,  in  spite  of  the  fact  that  it  is  not  quite 
a delta  fin.  (For  example,  compare  it  with  the  result  shown  in  Figure  39 
computed  by  slender  body  theory) . 

All  the  above  results  closely  agree  with  those  given  in  Reference  43, 
on  which  the  present  calculation  is  mainly  based. 

(d)  Cruciform  Slender  Delta  Fin-Body  Combination: 

Figure  55  shows  the  calculated  longitudinal  pressure  distribution 
on  the  body.  A pair  of  horizontal  fins  are  canted  by  3°  so  as  to  cause 
a negative  rotation  about  the  body  axis  (i.e.,  x-axis)  in  the  sense  of 
the  commonly  accepted  sign  convention.  It  should  be  noted  that  the  rate 
of  perturbation  pressure  change  is  very  sharp  at  the  leading  and  trailing 
edges  and  the  magnitude  of  the  pressure  is  lower  when  compared  with  the 
case  of  a large  aspect  ratio  fin  (see  Figure  52) , This  is  expected 
because  the  fin  is  slender  so  that  the  influence  of  the  fin  on  the  body 
should  be  small.  This  result  agrees  with  the  concept  of  slender-body 
theory  which  suggests  that  the  aerodynamic  disturbance  is  restricted  to 
the  cross  plane. 

4.  Discussions  on  the  Compressibility  Effect.  The  effect  of 
compressibility  on  the  lift  coefficient  of  a cruciform  slender  delta  fin- 
body  combination  is  shown  in  Figure  56.  (Only  the  angle  of  attack  influ- 
ence is  considered  here  for  the  sake  of  discussion) . The  fin  area  was 
extended  through  the  body,  as  shown  by  the  dotted  line,  in  order  to 
specify  a reference  area  for  the  lift  coefficient  computation.  In  general, 
the  Prandtl-Glauert  rule  results  in  an  over-correction  for  this  type  • 
small  aspect  ratio  fin  configuration  (see  Figure  56)  and,  [45], 

It  is  more  appropriate  to  use  Goethert's  rule  to  calculate  the  com- 
pressibility effect  on  fin-body  combinations.  For  a slender  fin-body 
combination,  after  using  the  Goethert  rule  correction,  the  result  (see 
Figure  56)  agreed  very  well  with  the  result  computed  by  the  slender 
body  theory.  In  view  of  a wider  range  of  Mach  number  applicability  of 
the  slender  body  theory,  such  an  agreement  is  not  totally  surprising..  In 
the  same  figure,  a result  due  to  Nielsen  [39]  is  also  included  for  com- 
parison. For  a non-slender  fin  case,  the  result  obtained  by  using  a 
Goethert's  rule  correction  is  shown  in  Figures  57  and  58. 
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Data. 


5.  turner's  Approach  Compared  with  Other  Methods  and  with 


(a)  Comparison  with  Slender-Body  Theory: 


A comparison  of  the  results  obtained  from  the  results  obtained  from 
the  singularity  method  and  from  slender-body  theory  has  previously  been 
briefly  discussed.  A surprisingly  good  agreement  between  the  two  methods 
was  obtained  for  the  lift  coefficient  of  a cruciform  slender  delta  fin- 
body  combination. 


The  pressure  and  load  distribution  on  the  same  fin-body  combination 
as  used  in  the  slender  body  theory  example  (see  Figure  25)  was  calculated. 
The  results  are  shown  in  Figures  59  to  61.  Figure  59  shows  the  pressure 
distribution  on  one  horizontal  fin.  It  is  worthwile  to  note  that  the 
linearized  expression  was  used  in  the  present  analysis  while  the 

expression  used  for  the  slender-body  theory  was  of  second  order.  The 
pressure  coefficient  was  computed  azimuthally,  at  six  positions  on  the 
one-half  body,  and  is  shown  in  Fieure  60.  Note  that  the  difference  in  the C 


expression  is  more  sensitive  on  the  fin  than  on  the  body.  This  means  that 
the  second  order  expression  for  the  is  more  important  for  such  slender 


fin-body  combinations  even  though  its  absolute  value  is  small.  Figure  61 
shows  the  load  distribution  on  both  a horizontal  and  a vertical  fin.  It 
can  be  seen  from  Figures  59  to  61  that  the  results  of  the  present  analysis 
and  of  slender-body  theory  agree  satisfactorily  with  each  other.  From 
these  comparisons,  it  can  be  concluded  that  both  methods  will  yield  a 
similar  result  for  a small  aspect  ratio  fin-body  combination.  This 
merely  proves  that  the  singularity  method  can  be  applied  for  fins  with  a 
wide  range  of  aspect  ratios  and,  therefore,  is  more  versatile. 


(b)  Compa'risons  with  the  "Flexstab"  Singularity  Distribu- 
tion Method  and  with  Data: 


The  "Flexstab"  method  [46]  has  been  applied  to  an. aircraft  wing-body 
combination  for  flight  at  high  subsonic,  transonic  and  supersonic  speeds 
[44  and  46],  A comparison  of  computations  based  on  K'orner's  scheme,  as 
adapted  in  this  analysis,  was  made  for  those  configurations  calculated 
by  the  "Flexstab"  technique. 

The  paneling  scheme  employed  for  the  present  calculations,  as  well 
as  those  by  the  "Flexstab"  method,  are  given  in  Figure  62.  The  main 
difference  between  the  "Flexstab"  method  and  Korner's  technique  is  in 
the  treatment  of  the  "body".  The  "Flexstab"  method  treats  the  body  by 
a number  of  discrete  vortices  distributed  on  the  mean  body  surface.  The 
mean  body  surface  is  approximated  by  a number  of  rectangular  strips  with 
streamwise  surfaces  [ 44] . In  addition,  the  "Flexstab"  method  is  more 
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complicated  and  takes  into  account  thickness,  camber  and  twist  effects 
on  the  wing  by  means  of  an  additional  source-sink  and  vortex  singularity 
distribution. 

The  calculated  chordwise  pressure  distribution  on  the  wing  is  shown 
in  Figure  63.  The  angle  of  incidence  of  the  body  was  taken  to  be,  2°,  and 
the  free-stream  Mach  number  was  0.85.  The  computation  was  made  near  the 
spanwise  mid-point.  The  agreement  with  "Flexstab"  results  was  surpris- 
ingly good.  The  agreement  with  the  experimental  data  were  also  satisfac- 
tory, except  in  the  region  very  close  to  the  wing  leading  edge,  where  a 
sudden  pressure  drop,  due  to  local  flow  separation,  is  found  in  the  data. 
The  case  for  an  angle  of  incidence  of  8°  is  show  in  Figure  64.  The 
agreement  with  the  "Flexstab"  method  was  very  good  in  this  case  also; 
however,  the  agreement  with  data  was  not  too  good  from  the  leading  edge 
to  the  mid-chord.  The  reason  for  this  disagreement  with  data  may  be 
attributed  to  flow  separation  from  the  sharp  leading  edge  because  of  the 
high  angle  of  incidence.  A further  study  for  such  a situation  is  needed. 

A further  comparison  with  the  "Flexstab"  method  was  made  for  the 
effect  of  a wing  flap.  The  calculated  chordwise  pressure  distribution, 
based  on  both  methods,  is  given  in  Figure  65.  In  this  case,  it  is  assumed 
that  the  wing  hasO0  angle  of  attack  and  that  the  flap  angle  is  8.3°.  The 
agreement  between  the  methods  is  good;  however,  the  inviscid  flow  theory 
predicts  a larger  pressure  difference  than  is  indicated  by  the  data  over 
the  flap  surface.  Again,  the  disagreement  is  probably  attributed  to  flow 
separation  on  the  upper  wing  resulting  from  the  high  flap  angle.  Finally, 
the  normal  force  coefficient  has  been  compared  as  shown  in  Figure  66. 

The  agreement  with  the  experiment  is  very  good.  The  compressibility 
effect  based  on  the  Goethert's  similarity  rule  is  again  proven  to  be 
very  satisfactory. 

D.  Comparisons  of  Fotential  Flow  Theories  with  Data  on  Slender 

Clipped  Delta  Fin-Body  with  Elevon  Deflection 

The  computational  model  employed  for  a cruciform  slender 
clipped  delta  fin-body  combination  is  given  in  Figure  67.  Each  fin  has 
a small  control  fin,  called  an  "elevon",  and  it  is  assumed  that  its  hinge 
line  is  located  near  the  trailing  edge  of  the  main  fin. 

The  angle  of  incidence  of  the  body  is  taken  to  be  8.54°  with  an 
elevon  deflection  angle  of  ± 30°,  so  that  the  body  generates  a counter- 
clockwiserolling moment  then facing  the  upstream  direction.  The  free- 
stream  Mach  number  is  1.62.  With  proper  sweptback,  as  in  this  example 
calculation,  the  leading  edge  of  the  fin  is  submerged  in  a subsonic 
stream.  Because  of  the  slenderness  of  the  fin  and  its  subsonic  leading 
edge,  one  expects  that  a better  agreement  would  result  by  employing 
slender  body  theory  than  from  use  of  the  singularity  distribution  method. 

A comparison  with  the  data,  indicates  that  this  is  indeed  true. 
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Figures  68  to  73  show  the  chordwise  pressure  distributions,  computed 
by  the  singularity  distribution  and  by  the  slender  body  theory,  at  dif- 
ferent spanwise  positions,  together  with  experimental  measurements  [45], 
The  agreement  of  the  singularity  distribution  method  of  KBrner  with  the 
data  was  adequate  in  general  over  the  fin  surface,  even  though  the  flow 
is  supersonic.  The  pressure  distribution  on  the  elevon  surface  could  not 
be  compared  due  to  a lack  of  data. 

As  was  expected,  the  computed  result  from  slender-body  theory  showed 
better  agreement  with  the  fin  data.  This  is  especially  so  on  the  lower 
surface  of  the  fin.  However  this  method  cannot  calculate  the  effect  of 
the  elevon  deflection  or  the  clipped  tip  effect.  'Harmans'  method  esti- 
mated a rather  large  pressure  distribution  on  the  fin  close  to  the 
hinge  line  of  the  elevon  as  a result  of  the  very  large  elevon  deflection 
angle.  The  large  pressure  increase  is  again  attributed  to  flow  separa- 
tion. This  can  be  seen  more  clearly  on  the  upper  surface  of  the  fin  near 
the  hinge  line.  The  tip  effect  caused  by  the.  leading  point  P of  the  tip 
Mach  cone  may  cause  the  pressure  to  increase.  However  no  correction  has 
been  made  for  this  in  the  present  analysis. 


The  longitudinal  C 


P 


distribution  at  a fixed  azimuth  angle  on  the 


body  has  been  computed  and  is  shown  in  Figure  74.  The  fin-body  combina- 
tion in  this  example  has  angle  of  incidence  of  7.4  and  a bank  angle  of 
2.8°.  The  C distribution  was  computed  from  a linear  summation  of  the 


C 


P 


resulting  from  separate  calculations  for  angle  of  attack  and  for  yaw 


angle.  The  results  are  as  shown  in  Figure  74.  (Also,  see  Figure  27  for 

results  for  the  individual  flows) . The  trends  of  the  C distribution  in 

P 


this  case  is  very  similar  to  the  one  discussed  previously  in  Figure  55. 


Various  compressibility  corrections  for  the  pressure  coefficient  at 
a sample  point  on  the  fin  are  shown  in  Figure  75.  An  angle  of  attack  of 
8.5  was  assumed  in  this  case.  Computed  results  using  the  Prandtl-Glauert 
and  the  Karman-Tsien  rules  showed  the  over-correction  discussed  earlier. 
The  total  lift  coefficient  in  this  particular  fin-body  combination  with 
zero  elevon  angle  is  shown  in  Figure  76.  The  maximum  body  cross-section 
area  was  taken  as  the  reference  area  for  the  lift  coefficient  calculation. 
The  contribution  of  the  body  to  the  total  lift  coefficient  is  very  small, 
and  may  be  neglected. 


III.  ANALYSIS  OF  VISCOUS  FLOW  AND  INTERACTION  WITH  PLUME 

The  basic  theoretical  statements  associated  with  the  determina' 
tion  of  the  viscpus  flow  are  made  in  the  following  section.  Two  basic 
flow  areas  need  treatment.  In  the  first  is  included  the  boundary  layer 
development,  both  laminar  and  turbulent,  over  the  body,  while  the  second 
relates  to  the  strong  interaction  in  the  base  flow  region. 
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So  far  the  development  has  been  restricted  to  a body  of  revolution 
when  the  cross- flow  component  in  the  boundary  layer  is  not  a feature. 

For  small  incidences,  when  it  is  reasonable  to  expect  only  minimal  cross 
flow,  then  it  is  possible  to  adopt  a strip  theory  approach  to  obtain  an 
approximate  solution. 

The  strong  interaction,  as  calculated  here,  is  essentially  a first 
order  result  since  no  attempt  has  been  made  to  correct  the  flow  over  the 
body  (ahead  of  the  interaction)  for  the  large  disturbance  created  by  the 
plume  interaction.  It  was  shown  earlier  [6]  that  in  the  absence  of  the 
plume  interaction  at  the  rear  of  the  body,  there  was  very  little  influence 
of  the  viscous  flow  on  the  potential  flow  structure.  While  with  a large 
plume  interaction  present,  it  is  likely  that  the  adjustment  to  the  fore- 
body flow  would  be  more  significant,  this  has  been  neglected  in  the  pre- 
sent study.  The  numerical  calculations  of  South  [47]  for  bodies  of 
revolution  show  the  expected  magnitude  of  this  correction.  In  a two- 
dimensional  flow  where  the  pressure  gradients  are  more  severe  and  the 
boundary  layer  growth  more  rapid,  it  is  shown  that  the  viscous  effect  is 
of  great  importance  [48], 

Figure  77  shows  the  overall  flow  field  for  the  type  of  situation 
under  discussion  and  identifies  the  principle  segments  of  the  flow.  A 
boundary  layer  develops  over  the  forward  part  of  the  body  (laminar  ini- 
tially followed  by  a turbulent  portion)  under  essentially  the  potential 
flow  pressure  distribution.  This  boundary- layer  may  pass  under  a shock 
interaction  and  may  negotiate  discontinuities  in  slope  at  a cylinder  boat- 
tail  junction  for  instance.  Near  the  rear  of  the  body,  the  boundary 
layer  development  is  influenced  by  the  strong  interaction  with  the  plume 
and  sufferes  both  a separation  process  and  a wake  confluence  before 
issuing  into  the  far  field. 

The  basic  concept  of-  the  methods  adopted  for  the  solution  of  the 
strong  interaction  shows  a close  connection  with  the  classical  component 
approach  [ 49] , but  there  have  been  certain  developments  and  modifications 
to  suit  the  situation  in  hand.  The  adoption  of  the  component  concepts  in 
no  way  suggests  its  superiority  over  an  integral  technique,  but  rather, 
was  chosen  on  the  basis  of  available  experience.  The  comments  made  later 
will  indicate  that,  in  the  present  writers'  opinion,  both  these  approaches 
have  been  developed  to  their  limit  and  are  still  not  completely  adequate 
in  all  situations.  There  is  a strong  urge  to  develop  numerical  methods 
for  the  entire  flow  field. 

Another  theme  runs  through  the  viscous  flow  development  as  presented 
here.  No  attempt  has  been  made  to  include  thermal  effects  in  the  analy- 
sis, which  implies  that  any  modification  of  the  strong  interaction  flow 
due  to  non-isoenergetic  considerations  is  neglected.  Base  heating  rates 
cannot  be  determined.  Consistent  with  this  overall  approach,  it  was  felt 
that  chemical  reactions  and  other  non-equilibrium  effects  could  be 
neglected  in  the  plume  calculations.  The  whole  theory  must  be  considered 
within  this  restriction. 
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A 


The  Viscous  Flow 


The  equations  of  motion  for  the  viscous  flow  over  a body 
of  revolution  can  be  stated: 


^ (Rpu)  + ^ (Rpv) 
du  , chi  1 c)p 


= o 


.2 

V-^-y-  + T 

c)r 


(43) 

(44) 


Over  the  leading  edge  of  the  body,  the  flow  is  laminar  and  then  the  con- 
tribution ' T 1 from  the  Reynolds  shear  stress  vanishes  identically.  R 
denotes  the  local  body  radius. 

The  viscous  flow  calculation  over  the  body  is  to  serve  two  major 
purposes:  determination  of  the  body  drag  and  evaluation  of  initial  con- 
ditions at  the  start  of  the  plume  induced  separation.  For  both  these 
activities,  it  is  only  necessary  to  determine  the  integral  characteristics 
of  the  boundary  layer  (the  momentum  thickness  and  skin  friction  being  the 
most  significant  quantities  in  this  connection) . 


With  this  in  mind,  it  was  considered  unnecessary  to  solve  the  boun- 
dary layer  equations  in  great  detail  and  large  numerical  expense.  Hence 
the  Thwaites'  method  was  adopted  for  the  laminar  flow  (in  conjunction  with 
Stewartson/Illingworth  and  the  Mangier  transformations)  and  Nash's  tech- 
nique was  used  for  the  turbulent  part  of  the  boundary  layer. 

Certain  limitations  implicit  in  this  approach  should  be  clearly 
recognized  at  the  outset.  Integral  methods  cannot  treat  flows  containing 
discontinuities  without  some  additional  empiricism.  Hence  the  boundary 
layer  development  over  such  configurations  as  a cone  or  boattail/cylinder 
junction  cannot  be  treated  without  special  attention.  The  same  is  true 
for  the  interaction  of  a shock-wave  with  the  boundary  layer  and  the 
individual  treatments  reserved  for  these  flows  will  be  dealt  with  in 
subsequent  sections. 

Throughout  the  boundary  layer  calculation,  the  external  pressure  has 
been  determined  from  the  potential  flow  solution  that  was  discussed  in 
Section  II  for  various  body  configurations.  .As  yet  there  has  been  no 
attempt  to  calculate  the  boundary  layer  development  over  the  fin  or  at 
the  fin-body  junction.  In  cases  where  plume-induced  separation  has  been 
calculated  for  finned  bodies  [24],  this  has  been  done  on  a strip  theory 
basis  with  those  segments  close  to  the  fin  root  being  neglected.  The 
wake  issuing  from  the  fin,  and  its  interaction  with  the  plume,  has  also 
been  disregarded,.  It  is  felt  that,  provided  the  fin  incidence  is  suf- 
ficiently low  to  assure  attached  flow  over  the  fin,  then  the  modifica- 
tion to  the  viscous  layer-plume  surface  mixing  will  be  minimal. 
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1.  The  Laminar  Boundary  Layer  Development.  The  theoret- 
ical treatment  of  the  laminar  boundary  layer  development  rests  upon  the 
application  of  two  transformations  to  the  calculation  procedure  of 
Thwaites  [50].  First,  there  is  the  Mangier  transformation  [51,  52]  which 
relates  the  flow  past  a body  of  revolution  to  that  past  an  equivalent 
two-dimensional  body,  then  the  Stewartson/Illingworth  transform  [53,  54] 
converts  the  result  to  a compressible  flow.  The  limitations  implicit  in 
the  use  of  these  transformations  (for  example,  the  boundary  layer  must 
be  thin  compared  to  the  body  diameter,  the  fluid  has  unit  Prandtl  number, 
etc.)  has  been  reviewed  elsewhere  [1  and  7]  and  need  not  be  restated 
here.  In  the  present  application,  these  restrictions  should  not  represent 
too  severe  a limitation. 

As  in  Reference  7,  we  can  state  the  final  result  for  the  laminar 
boundary  layer  momentum  thickness  development  in  the  form: 


(45) 


Here,  L is  the  characteristic  length  of  the  problem  (taken  to  be  the 
body  length  in  this  case)  and  R(x)  is  the  local  body  radius.  The  local 
velocities  Ue  and  the  local  temperature  Tg  are  simply  related  to  the 

local  Mach  number  when  the  external  flow  is  isentropic.  Tq  is  the  local 

total  temperature.  The  Reynolds  number  R^  used  in  this  equation  is  based 

upon  the  characteristic  length  L. 


The  solution  is  completed  by  employing  the  empirical  data  of  Thwaites 
[50]  to  determine  the  form  factor  H,  (and  hence  the  displacement  thick- 
ness) and  the  local  friction  coefficient,  C^. 


In  particular: 
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The  following  equations  were  fitted  to  Thwaites'  data: 

0.222  - 1.557  P - 1.781  P2 

= 0.222 

f 9 

0.222  - 0.466  P - 2.222  P 

•2.61  + 3.851  P + 5.7682  P2  + 0.4409  P3 
H = 2.61 


P < 0 
P = 0 
P > 0 


P < 0 
P = 0 


2.5728  + 14.2214  P - 432.396  P + 5132. 96P  P > 0 
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and  the  pressure  gradient  parameter  P is  given  as: 
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A typical  calculation  qf  the  laminar  boundary- layer  over  an  ogive 
forebody  is  shown  in  Figure  78.  Here  the  skin  friction  and  momentum 
thickness  variations  for  several  Reynolds  numbers  are  presented  for  the 
flow  over  an  ogive  forebody.  The  free-stream  Mach  number  is  1.1. 


2.  Transition.  No  attempt  has  been  made  to  predict  a 

transition  point  in  the  present  study.  Two  courses  are  open  in  the 
execution  of  the  calculation  procedure.  First,  a transition  point  can 
be  inserted  as  an  arbitrary  parameter  in  the  calculation  and  the  laminar 
boundary  layer  development  is  determined  up  to  that  point.  If  the  sepa- 
ration condition  (C^  = 0)  is  reached  during  the  progress  through  this 

laminar  flow  calculation,  then  the  calculation  is  stopped  and  transition 
is  assumed  to  occur  at  this  point.  The  remainder  of  the  boundary  layer 
is  determined  from  the  turbulent  boundary  layer  calculation  procedure. 

In  either  case,  the  final  value  of  momentum  thickness  at  the  termination 
of  the  laminar  flow  calculation  is  used  as  an  initial  condition  at  the 
start  of  the  turbulent  flow. 


3,  The  Turbulent  Boundary  Layer  Development.  The 
boundary  layer  equations  43  and  44  are  solved  by  first  integrating 
across  the  boundary  layer  to  give: 
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1 a / n \ dU 

- ~ (p  H 80  = t - p U S*  (48) 

R dx  \e  e / w Ke  e dx  N ' 

which  is  the  momentum  integral  equation  for  a turbulent  boundary  layer. 
This  equation  can  then  be  solved  numerically  provided  that  sufficient 
additional  information  concerning  the  physical  nature  of  the  turbulence 
is  included.  Nash  [55]  selected  various  empirical  data  to  close  the 
system  of  equations  for  a compressible  turbulent  boundary  layer. 

In  addition  to  Equation  48,  it  is  required  that  the  kinetic  energy 
integral  equatioh, 


_1_ 

2R 


d / 
& (p 


d / tt3 
[p  U 

k e e / 


W*  ) = 


Su 

by 


d y 


2 - e 

p u -r  &** 

e e dx 


is  also  solved. 


(49) 


The  variable  G is  introduced  to  determine  the  variation  of  the  tur- 
bulence production  integral.  Let 

G = i/2  C + — -j  Xc  T dy 

p u J 

Ke  e 

and  assume  that  the  function  G is  determined  from  the  differential 
equation. 

g -fco[/T  -/G]  (50) 

where  the  'constant'  G is  given  as 

G = -0.00883  + 0.02781  H - 0.03488  H2 

u u 

(5i) 

+ 0.01950  H3  - 0.00359  H4 
u u 


and,  the  "shape  factor"  Hu,  relates  to  a flow  with  no  density  changes 
(i.e.,  is  identical  to  H when  the  flow  is  incompressible). 

The  various  integral  quantities  can  be  related  if  the  velocity  pro- 
file relations 
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(52) 


are  adopted. 

Nash  [ 55]  showed  that  the  following  empiricism  is  adequate  for  the 
present  purposes: 
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This  completes  the  solution,  except  for  the  boundary  layer  thickness,  6, 
which  can  be  determined  from 


exp  [-3.5  (hu  - l)]  } 

1 + 0.1  M2 
e 


(54) 


The  differential  relation  expressed  by  Equation  50  goes  some  way 
towards  modeling  the  behavior  of  the  turbulence  structure.  Different 
values  of  the  constant  \ are  taken  as 


\ = 4 dG/dx  > 0 

X » 10  dG/dx  < 0 


so  that  different  boundary  layer  responses  are  incorporated  depending 
upon  the  nature  of  the  flow's  departure  from  equilibrium. 
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In  simple  physical  terms,  this  differential  equation  just  assesses 
the  rate  at  which  the  turbulent  boundary  layer  would  return  to  the 
equilibrium  condition  after  suffering  some  perturbation.  While  this  is 
clearly  superior  to  a simple  algebraic  relation  for  departure  from  equili- 
brium, it  still  does  not  incorporate  very  much  detail  about  the  physical 
structure  of  the  Reynolds  stress  tensor  and  its  response  to  an  external 
distrubance.  In  consistency  with  the  basic  premise  that  only  an  indica- 
tion of  the  integral  properties  of  the  boundary  layer  are  required  in  the 
present  application,  it  was  considered  that  such  a simple  representation 
of  the  turbulent  stresses  would  be  adequate. 

The  turbulent  boundary  layer  calculation  is  started  at  the  transition 
point  with  an  initial  value  of  momentum  thickness,  0,  equal  to  that  at 
the  termination  of  the  laminar  flow  calculations.  The  initial  value  of 
H was  taken  to  be 

H - 1.6  + 0.417  M^ 
et 

where  Mefc  is  the  local  Mach  number  at  the  transition  point. 

Calculated  skin  friction  distributions  are  displayed  on  Figure  79 
for  three  different  Reynolds  numbers.  The  curves  show  the  change  in 
skin  friction  at  the  transition  point  and  also  at  the  boattail  shoulder. 
Details  of  the  calculation  procedure  at  the  boattail  slope  discontinuity 
are  discussed  in  Subsection  4 of  the  current  section. 

4.  Change  in  Integral  Parameters  Across  Discontinuities. 
Several  flow  situations  of  interest  require  knowledge  of  the  turbulent 
boundary  layer  development  across  some  form  of  discontinuity.  Particu- 
larly, this  involves  the  interaction  of  a weak  shock  wave  with  a boundary 
layer  under  conditions  that  do  not  involve  separation,  and  the  flow  over 
a slope  discontinuity.  The  latter  situation  develops  at  the  junction  of 
a boattail  or  flare  or  at  the  shoulder  of  a conical  nose. 

In  the  current  development,  it  is  considered  that  only  a very  appro- 
ximate analysis  is  appropriate  since  the  entire  boundary  layer  calculation 
is  attempted  within  the  context  of  the  integral  theory.  Not  even  the 
more  sophisticated  approaches  to  such  flows  [56]  are  adopted,  although 
they  could  be  with  little  effort. 

Basic  to  the  approach  adopted  is  the  idea  that,  if  the  pressure 
change  is  sufficiently  rapid,  then  any  changes  that  it  may  produce  in 
the  turbulence  structure  can  be  neglected  to  the  first  order  and  changes 
in  the  integral  parameters  determined  from  an  inviscid  model.  The 
experimental  work  of  Rose  and  Johnson  [57]  for  the  interaction  of  an 
oblique  shock-wave  with  a turbulent  boundary  layer  tends  to  confirm  this 
flow  model:  in  the  outer  (wake)  region  of  the  boundary- layer  the  turbulence 
structure  remains  more  or  less  frozen.  As  is  often  observed  in  such 
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interactions  in  turbulent  flow  158],  however,  the  wall  region  quickly 
returns  to  equilibrium.  This  is  to  be  expected  when  the  pressure  gradient 
across  the  viscous  region  is  small  the  wall  boundary  conditions  must  be 
imposed. 

Consistent  with  the  basic  model  (and  the  further  requirement  that 
the  total  change  be  small) , Nash  [ 59]  developed  the  relation: 

p U H 0 = constant  (55) 

Ke  e e 

as  the  first-order  change  in  momentum  thickness  across  a discontinuity. 

In  the  present  work,  this  relationship  is  adopted  in  the  treatment  of 
both  expansive  and  compressive  pressure  discontinuities. 

The  form  factor,  H,  of  the  boundary  layer  was  also  allowed  to  change 
discontinuously.  The  magnitude  of  the  change  being  given  by: 
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where  the  suffix  2 denotes  downstream  conditions  at  the  boundary  layer 
edge. 


New  values  of  0 and  H thus  allow  the  turbulent  boundary  layer  calcu- 
lation to  continue  and  other  quantities  (friction  coefficient  say)  are 
returned.  The  complete  boundary- layer  development  up  to  the  start  of 
the  plume  interaction  is  now  available. 


In  some  situations,  particularly  where  the  flow  downstream  of  the 
discontinuity  is  subject  to  an  adverse  pressure  gradient,  the  subsequent 
flow  development  can  be  sensitive  to  initial  values  as  calculated  from 
the  above  relations.  Figure  80  shows  this  effect  for  the  boundary  layer 
calculated  on  a conical  boattail.  The  potential  flow  theory  (Section  II) 
for  this  geometry  yields  a discontinuous  decrease  in  pressure  at  the  boat- 
tail  junction.  In  reality  there  would  be  a smooth,  but  rapid,  change. 

The  calculations  reported  on  Figure  80  were  for  different  estimates  of 
this  smoothing  of  the  pressure  change  across  the  surface  slope 
discontinuity. 


Based  on  this  result  it  appeared  that  a further  study  of  the  problem 
was  desirable  and  an  experimental  program  was  set  up  to  study  the  flow 
development  over  slope  discontinuities.  A typical  two-dimensional  model 
is  shown  mounted  in  the  wind  tunnel  in  Figure  81.  More  details  of  the 
wind  tunnel  are  given  in  Reference  60.  The  slope  discontinuity  afforded 
by  this  model  configuration  was  6,5°  and  its  length  was  8.8  inches. 
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Surface  static  pressures  were  measured  over  the  entire  model  length 
and  boundary  layer  velocity  profiles  were  obtained  at  specific  locations. 
Figure  82  shows  a measured  surface  static  pressure  distribution  at  a 
free- stream  Mach  number  of  0.9.  As  expected,  the  flow  suffers  a strong 
expansion  at  the  corner  followed  by  a recompression.  Due  to  the  boundary 
layer  development  over  the  configuration,  the  flow  does  not  recompress 
to  a stagnation  condition.  Far  downstream  it  is  found  that  the  flow 
returns  to  the  free- stream  pressure. 

Typical  mean  velocity  profiles  are  presented  on  Figure  83.  It  is 
shown  that,  as  the  flow  approaches  the  expansion  corner  from  the  undis- 
turbed upstream  (flat  plate)  boundary  layer,  there  is  a marked  thinning 
of  the  boundary  layer  with  a consequent  increase  in  wall  shear.  After 
the  flow  has  negotiated  the  rapid  expansion  at  the  corner,  the  boundary 
layer  thickens  and  the  wall  shear  is  reduced. 

The  integral  techniques  for  calculating  the  turbulent  boundary  layer 
development  do  not  include  details  of  this  type  when  applied  to  flow  over 
slope  discontinuities.  Neither  the  calculation  techniques,  nor  the  experi 
mental  study  considered  here,  discussed  the  change  in  turbulence  structure 
that  is  associated  with  the  flow.  This  subject  must  receive  greater 
attention  in  the  future. 


A somewhat  different  problem  exists  at  the  junction  of  a flared  rear 
end  with  the  cylindrical  forebody.  The  potential  flow  theory  discussed 
in  Section  II  would  place  a stagnation  point  at  this  junction  and  does 
not  in  any  way  reflect  the  real  flow  development.  If  this  potential  flow 
pressure  distribution  were  adopted  as  the  starting  point  for  a boundary 
layer  calculation,  then  it  is  evident  that  a separation  would  be  predicted 
and  the  calculation  terminated.  The  experimental  data  of  Reference  61. 
was  evaluated  to  determine’  an  empirical  relation  for  this  junction  pres- 
sure so  that  the  boundary  layer  calculation  could  continue  uninterrupted 
to  the  end  of  the  body.  With  the  notation 
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where  P^/P^  is  the  static  pressure  ratio  at  the  flare  junction  and  0 ^ 

the  flare  angle,  the  correlation  is  as  shown  in  Figure  84,  and  is 
included  in  the  computer  program.  Details  of  the  body  geometries  used 
in  the  studies  are  summarized  in  Appendix  B. 


5.  Calculation  of  the  Free-Shear  Layer.  Several  restric 
tions  from  complete  generality  are  contained  in  the  free- shear  layer  cal- 
culation. First,  it  is  assumed  that  the  flow  is  turbulent  so  that 
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no  comparison  between  laminar  and  turbulent  base  flow  is  possible.  It 
is  also  assumed  that  the  shear  layers  evolve  along  a conical  trajectory 
and  that  their  profile  can  be  approximated  by  the  asymptotic  turbulent 
shear-layer  profile.  Two  major  assumptions  are  involved  here.  First, 
it  is  supposed  that  the  effect  of  the  initial  velocity  profile  (i.e.,  that 
associated  with  the  separating  boundary  layer  at  the  start  of  the  inter- 
action) is  insignificant  in  as  far  as  the  distribution  of  mass  flow  in  the 
shear  layer  is  concerned.  This  latter  quantity  is  significant  since  the 
solution  to  the  base  flow  problem  is  closed  by  means  of  a mass  balance 
in  the  dead  air  region.  An  approximate  correction  for  this  initial  pro- 
file effect  is  included  [ 61]  by  definition  of  a coefficient  of  the  form: 


Coeff  = 1.21  - 0.004  x/0  x/0  < 50 


= 1.00  ' x/0  > 50 


(58) 


which  changes  the  location  of  the  median  streamline  (see  Equation  66) . 
The  second  assumption  relates  to  the  neglect  of  any  external  influences 
on  the  shear  layer  development.  In  this  connection  it  is  supposed  that 
the  external  stream  is  at  constant  pressure  and  that  the  shear  layer 
development  is  not  influenced  by  its  close  proximity  to  the  body  surface 
or  to  the  jet  plume.  It  is  certainly  not  true  that  the  external  flow  is 
a constant  pressure  flow,  but  the  relative  shift  between  the  streamline 
dispositions  resulting  from  the  pressure  gradient  and  the  upstream  pro- 
file effect  is  not  known. 


Although  the  shear  layer  velocity  profile  does  not  match  that  of  the 
approaching  boundary  layer,  this  is  acknowledged  to  the  extent  that  the 
momentum  thickness  of  the  shear  layer  is  selected  to  be  equal  to  that  of 
the  boundary  layer.  Two  cases  must  be  considered.  When  the  base  pressure 
is  lower  than  the  ambient  pressure  then  the  boundary  layer  must  negotiate 
an  expansion  fan  at  the  base  of  the  body.  The  equations  discussed  in 
Subsection  3 provide  the  change  in  boundary  layer  momentum  thickness: 

^e^e^e®  = Constant  (59) 


so  that  the  value  at  the  start  of  the  free-shear  layer  can  be  estimated. 
This  relation  is  also  used  at  the  nozzle  exit  to  determine  the  change  in 
boundary  layer  momentum  thickness  at  the  start  of  the  plume  shear  layer. 
The  other  case  occurs  when  the  boundary  layer  leaves  the  body  ahead  of 
the  base  plane  due  to  the  interaction  of  a shock  wave.  Now  the  analysis 
of  Stratford  [62]  and  McDonald  [63]  shows  that 


1 + 


7-1  m2 


1 3 


0SL  = 6BL 


1 + ^T1 


(60) 
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where 


A - 1.0  - 0.605 


Here,  M0  denotes  the  approach  Mach  number  at  the  start  of  the  interaction 

and  denotes  the  Mach  number  downstream  of  the  shock  wave.  This  latter 

value  is  determined  from  the  theory  of  Mager  [64]  and  is  discussed  in 
Subsection  1 

Finally,  it  is  assumed  that  the  entire  shear  layer  development  takes 
place  under  isoenergetic  conditions  and  base  heating  effects  and  hot 
exhaust  plumes  are  not  treated. 

With  these  assumptions,  it  is  well  known  that  the  velocity  profile 
is  given  in  the  form: 

u/Ue  = 1/2  £l  + erf  (t])J  (61) 

where  the  similarity  parameter  p is  p = o y/x.  The  coordinate  among  the 
shear  layer  is  x and  y is  perpendicular.  The  mixing  coefficient,  a,  will, 
be  related  to  the  Mach  number  at  the  edge  of  the  shear  layer  (M  ) by  the 
Korst  and  Tripp  [ 64]  relation 

a = 12.0  (1  + 0.23  M ) (62) 


As  it  stands,  this  error  function  velocity  profile  is  for  a two- 
dimensional  shear  layer  but  can,  with  restriction,  be  adopted  for  the 
axisymmetric  flow  situation.  Provided  that  the  boundary  layer  remains 
thin  compared  with  the  diameter  of  the  body  and  that  the  confluence 
point  is  not  located  too  close  to  the  axis  of  symmetry  then  the  Mangler- 
Stepanov  transformation  could  be  adopted  [51,  52]  . It  should  be  cau- 
tioned, however,  that  this  transformation  is  only  valid  for  a laminar 
flow  and  its  application  to  the  turbulent  shear  layer  is,  at  best,  an 
approximation  [ 65]  . 


In  the  present  implementation  of  the  transformation  it  is  supposed 
that  the  shear  layer  is  essentially  conical  so  that  the  modified 
similarity  parameter  is 


3R  R a 

o p 

n = 2 2 

R R + R + R 
op  o p 


X 


(63) 
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where  R is  the  shear  surface  radius  at  the  base  plane,  R is  the  value  at  the 
o P 

station  of  interest  and  x and  y are  the  coordinates  in  the  axisymmetric.  con' 

a a 

figuration.  The  value  of  the  mixing  coefficient,  a,  is  left  unchanged 
without  justification. 

The  streamlines  within  the  shear  layer  are  located  on  the  basis  of 
Nash's  development  [59],  Then  it  is  found  that; 


(64) 


where  x is  the  coordinate  along  the  streamline  and  0 is  the  value  of 
momentum  thickness  [from  either  Equation  (59)  or  Equation  (60)  and  the 
appropriate  boundary  layer  calculation]  at  the  start  of  the  shear  layer. 
As  shown  in  Reference  61,  the  function  G is  expressed  by  the  equation: 


G,<n)  - 


1.0 


/ 


exp  (q  ) (u/Ue) 


a 


1 + zrMJ 


t1-  t] 


d(u/Ue) 


(65) 


with  the  velocity  profile  being  given  by  Equation  61.  The  velocity, 

Um/Ue>  along  the  median  streamline  was  given  [59]  to  be 

u _ -i  1/2 

j2  = 0.348  + 0.018  Ml 

and  determine  the  reference  function  G...  in  terms  of  the  external  Mach 

M 

number  M£.  Typical  distributions  of  the  function  G across  the  shear 
layer  are  shown  in  Figure  85. 

6.  Treatment  of  the  Strong  Interaction.  Free  shear 
layers  are  placed  along  the  surface  of  the  plume  and  as  a boundary  between 
the  external  flow  and  the  dead-air  region  (see  Figure  77).  The  solution 
is  closed  at  the  confluence  point  which  is  so  selected  as  to  give  a net 
mass  entrainment  in  the  dead-air  region  equal  to  zero. 

Equation  (64)  provided  a relationship  between  the  length  of  the 
shear  layer  and  the  initial  momentum  thickness,  with  the  external  Mach 
number  (assumed  constant  along  the  entire  shear  layer)  acting  as  a 
parameter.  The  two  shear  layers  are  matched  via  the  geometric  constraints 
of  the  problem  and  a pressure  criterion  at  the  confluence  point.  Now 
Equation  64  gives, 
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(67) 


and  the  condition 

Ats  A^s  + Af  = 0 


must  be  imposed  on  the  solution  to  satisfy  the  zero  net  entrainment  require- 
ment. In  these  equations,  the  suffix  " s " refers  to  the  free-shear  layer 
and  "p  "refers  to  the  shear  layer  along  the  plume  boundary.  The  momentum 
thickness  at  the  start  of  the  shear  layer  is  again  denoted  by  Q. 

Figure  86  shows  that  the  following  geometric  constraints  must  hold: 


y + R - R 

tan  (0  - 0 ) = IL 

s 7 1-x+x+x 

s n p 

2 2 2 
and  £ = (1  - x + x + x)  +(y  + R R ) 
s n p '7p  n s 


(68) 


2 2 2 
S = y + x 
P P 


provided  that  it  is  assumed  that  the  plume  boundary  is  conical  between 
the  nozzle  exit  plane  and  the  confluence  point. 


The  final  pressure,  Pc,  at  the  end  of  the  confluence  process  is 

unknown.  It  is  supposed  that  the  individual  pressure  rises  for  either 
the  external  or  the  internal  flows  can  be  expressed  in  the  form  of  the 
small  turning  relation  modified  by  an  undetermined  factor.  Thus  it  is 
supposed  that: 


P 

p2  - 1 + £(MS)  Tas  (69) 

s 


where  the  coefficient,  r,  is  determined  empirically,  and  relates  to  the 
N factor  of  Reference  61.  The  function  f(M)  has  the  form: 


f(M)  = 


yM 


yv 
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and  is  applied  for  either  subsonic  or  supersonic  Mach  numbers.  The  angle 

a (and  the  corresponding  angle  a ) are  defined  in  Figure  87.  The  local 
s p 

plume  slope  0^  is  determined  from  the  plume  geometry  and  then  ag  and 
are  related  by  the  equation: 

as  =ap  + 0L  -'0s  + 0-  (70) 


The  velocity  um  is  determined  from  Equation  66  and  hence  the  function 

Gm  is  known.  Now  the  only  remaining  unknown  quantity  is  the  velocity  on 

the  reattaching  streamline  which  determines  G via  Equation  (65) . It  is 
* ” 
assumed  that  the  flow  on  the  reattaching  streamline  is  reduced  to  a stag- 
nation condition  at  the- confluence  point  where  the  pressure  is  Pc>  The 

fact  that  this  pressure  is  unknown  accounts  for  the  quantity  r in  the 
relation  between  the  pressure  ratio  and  turning  angle  contained  in 
Equation  (69) . 


If  c = 


2+1  2 
2 M 


11/2 


H + 


zzi  M2 


Then  for  isentropic  compression  along  the  reattachment  streamline  to  a 
pressure  ratio  of  pc/pref  demands  that: 


(71) 


The  solution  procedure  is  as  follows.  An  angle  0 is  selected  (or 

s 

when  the  base  pressure  is  equal  to  the  plateau  pressure,  the  separation 
point,  xg,  is  the  parameter  since  then  0g  = and  the  above  equations 

solved  iteratively  to  determine  At«  The  initial  conditions  (0  or  x ) 

* s s 

are  systematically  changed  over  a range  of  values  until  Aty  changes  sign* 
Interpolation  then  gives  the  solution,  A^  = 0,  and  all  parameters  of  the 
base  flow  are  determined.  Of  particular  interest  are  the  base  pressure 
and  the  length  of  the  separation  region. 


Comparison  of  the  theoretical  calculation  with  experimental  data 
(a  typical  variation  of  base  pressure  with  r is  shown  on  Figure  88) 
indicated  that  it  was  appropriate  to  take 
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r = 1.277739  X M -0.356621  X M2  -0.69188:  z > 0.8 
s s ’ 

r = 0.1  z < 0.8 

7.  The  Separation  and  Plateau  Pressures.  When  the  base 
pressure  is  lower  than  the  ambient  pressure  it  is  possible  to  introduce 
an  expansion  fan  at  the  base  of  the  body  to  determine  the  conditions  along 
the  outer  edge  of  the  free-shear  layer  and  complete  the  calculation  using 
the  analysis  of  the  previous  section.  A different  situation  exists  when 
the  base  pressure  is  greater  than  the  ambient  value.  Now  it  becomes 
necessary  to  estimate  the  pressure  rise  that  occurs  under  the  compression 
in  the  external  flow. 

The  method  adopted  is  as  follows.  As  the  base  pressure  increases 

from  the  free- stream  value,  it  is  assumed  that  the  external  flow  deflects 

through  an  angle  0 such  that 
s 

9 y M2 

Pb/P  = 1 + (72) 


where  M is  the  local  Mach  number  at  the  start  of  the  interaction  (for  a 
parallel  sided  body  this  will  usually  coincide  with  the  free- stream  Mach 
number) . This  compression  process  cannot  continue  when  the  angle  0 

s 

becomes  too  large  since  then  the  boundary  layer  will  suffer  a separation 
from  the  body  surface.  In  the  present  analysis  it  is  assumed  that  the 
limit  point  at  which  flow  -separation  is  initiated  can  be  described  by  the 
theory  of  Mager  [34] . Details  of  this  theory  were  given  in  Reference  1 
and  need  not  be  repeated  herein.  It  may  be  remembered,  however,  that  this 
theory  does  not  include  any  reference  to  the  Reynolds  number  in  specifica- 
tion of  the  separation  pressure.  For  completeness,  we  note  that  the 
separation  pressure  [34]  is  given  by  the  relation: 


P /P„ 
s o 


= 1 + 


Zli 

2 


M2  K 
o 


1 + 


(73) 


where  the  constant  K has  been  given  the  value  0.45  as  the  best  representa- 
tion of  the  available  data.  Details  of  the  plateau  pressure  (which  pro- 
vides the  maximum  value  0^  to  the  flow  separation  angle  0g)  are  presented 

along  with  the  discussion  of  Equation  76.  The  Mach  Number  downstream  of 
the  interaction,  M2,  is  related  to  the  approach  Mach  number  MQ  by: 

M2  - /K  M0 

and  the  momentum  thickness  change  (Equation  60)  can  be  determined. 
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As  the  base  pressure  increased  from  the  ambient  value,  it  is  assumed 
that  the  flow  still  leaves  the  body  at  the  base  plane  until  the  base  pres- 
sure has  risen  to  the  separation  pressure  level.  At  this  condition  the 
separation  point  is  ready  to  move  forward  on  the  body,  but  it  is  still 
possible  for  the  base  pressure  to  increase.  It  will  continue  to  do  so 
until  the  separation  has  migrated  sufficiently  far  from  the  base  for  the 
entire  base  area  to  be  influenced  by  the  plateau  pressure. 


To  complete  the  solution  in  this  phase  of  the  flow  development,  it 
is  necessary  to  specify  both  the  plateau  pressure  and  the  streamwise 
length  scale  for  the  pressure  rise  as  a whole.  It  is  assumed  here  that 
the  pressure  signature  under  the  interaction  is  unique  and  independent 
of  the  position  of  the  separation  point  relative  to  the  base.  With  this 
assumption,  it  is  possible  to  adopt  a correlation,  somewhat  similar  to 
that  suggested  by  Zukoski  [ 66]  for  the  pressure  rise  through  the  inter- 
action. The  form  adopted  is  shown  in  Figure '89  and  the  two  length  scales 
xp  and  xq  required  in  this  correlation  are  shown  in  Figure  90.  This 

latter  curve  reduces  to  Zukoski' s values  when  the  approach  Mach  number  is 
greater  than  1.5  but  contains  an  elongation  of  the  interaction  length  at 
lower  Mach  numbers  roughly  in  accord  with  the  measured  distributions  of 
Rubin  [ 60] . The  extended  separation  length  is  also  used  for  subsonic 
approach  conditions  as  will  be  discussed  later. 


The  pressure  signature  of  Figure  90  is  fitted  by  the  relations 

Pp/Po  = 1,0  + (VPo  " 1,0)X  (74) 


where  the  function  X is  expressed  as 

X = 5 z2;  0 < z < 0.2 

X =»  Sin  [ it  (z  - 0.1)/1.8]  ; 1.0  > z > 0.2 

X = 1.0  ; z > 1.0 


(75) 


Here  z = (x  - xg  + xQ)/xp  and  xq  represent  the  location  of  the  separation 
point  and  the  length  scales  xg  and  x are  given  in  Figure  90. 


The  above  equations  specify  the  entire  pressure  rise  from  the 
approach  value  Pq  to  the  plateau  value  Pp.  Some  experimental  data  for 

the  plateau  pressure  are  shown  in  Figure  91  along  with  the  theoretical 
values  of  Mager  [ 34]  and  Zukoski  [ 66] . In  the  present  work  the  curve 
labeled  (A)  given  by 

Pp/PQ  = 0.915  + 0.484  Mq  (76) 
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has  been  adopted  for  specifying  the  plateau  pressure.  The  limiting  flow 
turning  angle  is  now  determined  from  Equation  72. 


Figure  92  indicates  how  the  separation  and  plateau  conditions  are 
situated  in  terms  of  the  development  of  the  base  pressure  signature  as 
a function  of  jet  pressure  ratio.  The  limiting  base  pressure  at  large 
jet  thrust  is  also  defined.  It  should  be  noted  that  much  of  the  total 
pressure  rise  above  the  ambient  value  takes  place  before  the  separation 
point  encroaches  upon  the  body. 

The  correlation  shoxm  by  Equation  (76)  is  only  adequate  for  two- 
dimensional  flows  and  for  flow  over  parallel-sided  bodies  of  revolution. 
When  the  rear  end  of  the  body  has  a conical  configuration  then  the  data 
of  Reference  62  shows  that  the  limiting  base  pressure  can  be  a strong 
function  of  conical  angle.  Only  limited  data  are  available  since  in  most 
of  the  cases  contained  in  Reference  62,  the  jet  pressure  ratio  was  not 
increased  sufficiently  to  reach  the  limiting  base  pressure.  Some  data 
are  presented  on  Figure  93  and  shows  the  trends  exhibited.  Similar  data 
for  two-dimensional  supersonic  boattail flow  are  given  in  Reference  67. 


A simple  idea  evolves  the  from  consideration  of  the  data  shown  in 
Figure  93.  It  is  suggested  that  the  ratio  of  the  limiting  base  pressure 
at  a given  conical  angle  to  that  at  0 angle  (P  /P  say)  is  almost  indepen- 
dent of  Mach  number.  This  reduction  of  the  data  is  shown  in  Figure  94 
and  substantiates  the  suggestion  for  supersonic  flow,  when  the  relation 


P /P  = 1 - 0.02440  8 + 0.00074  e‘ 
P 


correlates  the  data  well I All  that  remains  now  is  to  scale  the  plateau 
and  separation  pressure  relations  contained  in  Equations  73  and  76  to 
yield  more  realistic  limiting  conditions  for  conical  aft  ends. 

Finally,  it  is  of  interest  to  compare  the  predicted  pressure  rise 
under  the  separation  (from  Equation  74  and  the  related  structure)  with 
experimental  data. 

The  curves  on  Figure  95  show  the  theoretical  results  for  jet  pressure 
ratios  of  87.7  and  100  compared  with  the  experimental  data  [62] 
for  a jet  pressure  ratio  of  87.7.  The  free- stream  Mach  number  is  1.1. 

It  is  clear  that  the  theoretical  pressure  signature  is  of  the  correct 
form  and  the  length  scale  is  more  or  less  of  the  correct  order. 


When  the  approach  Mach  number  is  less  than  sonic  then  the  situation 
is  not  so  easily  predicted.  The  experimental  data  shown  on  Figure  93(b) 
suggests  that  the  separation  pressure  is  given  as, 
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P /P  = 0.565M  +0.795;  0.5  < M <1.0 

s co  o’  o 


The  corresponding  maximum  pressure  is  given  as: 

P /P  = P /P  + 0.01 

P CO  s oo 

and  the  difference  between  the  separation  and  "plateau"  conditions  is 
very' small. 


8.  Evaluation  of  Force  Coefficients.  Two  basic  drag 
integrals  are  required.  One  of  these  produces  the  form  drag  while  the 
other  generates  the  skin  friction.  For  the  purposes  of  analysis  it  was 
expedient  to  evaluate  the  form  drag  separately  for  each  component  of  the 
body. 

The  expression  for  the  base  drag  is  given  as: 

c = (pc° : pb)  ab 

DB  1/ 2p  U2  A 

roo  oo  ref 

■ 7 (1-VP-) 

In  the  present  case  the  reference  area  for  the  drag  coefficient  was  taken 
to  be  that  of  the  parallel  section  of  the  body. 

determined  from  the  equation 
dx  R(x)d0 


from  the  expression 
dx  R(x)d0 


friction  drag  the  approximation  that  the  skin 
friction  is  0 downstream  of  the  separation  point  has  been  made  and  the 
integration  is  stopped  at  the  separation  point.  The  form  drag  has  been 
broken  into  two  parts  (in  addition  to  the  base  drag) , and  the  boattail 
contribution  is  displayed  separately. 


The  friction  drag  is 


and  the  form  drag  follows 


In  the  calculation  of  the 


2 2 
h ' h 

M2  R2 
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The  usual  expression  for  a form  drag  involves  an  integral  with 
respect  to  the  body  radius: 

[ C dR 

J P 

In  situations  where  both  the  body  radius  and  the  pressure  are  varying 
rapidly,  it  is  more  accurate  to  adopt  the  following  as  the  appropriate 
integral: 

fc  ? dx 

J P dx 

provided  that  the  body  geometry  is  sufficiently  well  known, 

coefficient  are  expressed  by  the 


dx  R(x)  d 0 


xdx  R(x)  d 9 


nose  up  and  is  taken  about  the  nose 
of  the  body.  These  integrals  are  evaluated  on  the  basis  of  a strip  theory 
using  7 asithumal  locations. 

9,  Summary  of  the  Base  Flow  Calculation.  The  two  prin- 
cipal modes  of  base  flow  that  have  been  discussed  correspond  to  base  pres- 
sures either  below  or  above  the  ambient  pressure.  Figures  96  and  97  show 
these  two  flow  situations  and  summarize  the  equations'  required  for  their 
solution.  The  computer  program  (see  Appendix  A)  resolves  the  differences 
between  the  two  flows  without  operator  intervention. 

B.  The  Plume  Calculation 

If  it  is  assumed  that  the  jet  expands  inviscidly,  then 
the  plume  boundary  can  be  calculated  quite  accurately  by  the  method  of 
characteristics.  A schematic  graph  of  the  interaction  is  given  in 
Figure  98.  A short  development  of  the  characteristic  equations  used  for 
the  computation  of  the  plume  boundary,  assuming  that  the  pressure  distri- 
bution along  the  plume  varies  in  a given  way  follows  below: 

Consider  an  axisymmetric  flow  field  whose  axis  of  symmetry  coincides 
with  the  x-axis  of  a cylindrical  coordinate  system  (x,r,0)  as  shown  in 


The  normal  force  and  momenl 
equations: 


CN  = 


°M  “ 


* L II 

t^tJI 


C cos  9 
P 


C cos  0 
P 


where  the  pitching  moment  is  poi 
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Figure  99.  An  axi symmetric,  supersonic,  and  inviscid  flow  is  fully 
described  by  the  continuity  equation: 


1 

P 


Dt 


+ V.t  = 0 
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the  equation  for  the  conservation  of  momentum 


P 


_ _ 
Dt 


vp 


J 


(77) 


(78) 


the  rate  of  change  of  entropy,  s, 


(79) 


and  the  equation  of  state,  which  can  be  expressed  as 

P = p(P,s)  , (80) 

— > 

where  q is  the  velocity  vector,  p the  pressure,  and  p the  density. 
Equation  (79)  is  valid  everywhere  except  for  passage  through  a shock  waVei 

— > 

Equations  (77)  through  (80)  are  4 equations  for  the  4 unknowns  q, 
p,  p,  and  s and  can  be  combined  to  give: 


Ou  Ov  \ 
c)r  + Ox  J 


(81) 


written  for  an  axisymmetric  system  (x,r).  Now  consider  a line  in  the 
x,r,  plane  passing  through  a point,  P,  of  the  flow  inclined  to  the  x 

axis  by  the  angle  tan  then  the  expression  for  this  line,  which  is 
usually  called  a characteristic,  is 


(82) 


The  variations  of  the  velocity  components,  u and  v,  along  this  line  are 


du  _ Ou  Ou 
dx  Ox  Or 


(83a) 


dv  Ov  Qv  , 

dx  Ox  Or 


(83b) 


Substitution  of  Equation  (83)  into  Equation  (81)  and  using  the  expression 
for  irrotational  flow, 
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(84) 


3v  Su 
Sx  Sr 


gives 


+ (2K-HN.)  “ + (HX2  - 2KX+L)|j  + 7=0 


where 
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(85) 
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Now  consider  lines  through  the  point,  P,  for  which  the  equation 


HX  - 2KX  + L = 0 


(87) 


is  valid.  This  quadratic  equation  for  the  tangent  has  2 solutions  given 
in  the  form: 


and  the  2 lines  through  the  point  P are  defined  by 
dr 


dx  1 


dx  ^2 


Equation  (85)  becomes 


iH.v  iv  1 v = 

dx  2 dx  Hr 


(88) 


(89) 


(90a) 


and 


ia+x  f£  + i2. 0 

dx  1 dx  Hr 


(90b) 


The  flow  directions  defined  by  Equation  (88?are  real  only  if  the  inequality 

(91) 


K2-HL=M2-1>0 


is  satisfied  and  the  local  velocity  is  greater  than  that  of  sound. 

The  expressions  for  the  tangents  given  in  Equation  88  can  be 
simplified  to  read 


Xl  2 = tan  (9  - l1) 


(92) 
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where  0 is  the  inclination  of  the  velocity  vector  and  p.  is  the  local 
Mach  angle.  Hence,  the  expressions  (90a)  and  (90b)  are  valid  along  the 
lines  whose  tangents  are  inclined  at  the  Mach  angle,  p , with  respect  to 
the  local  streamline  direction.  These  are  the  characteristic  lines  for 
the  flow.  Similar  to  the  treatment  of  Equation  (92),  Equations  (90a)  and 
(90b)  can  be  simplified  by  the  introduction  of  the  variables  p and  0 to 
give: 


cot  p 


ia_de 

q 


smu  sin6  , 

7z — ; — ■ dx 

cos(0+p)r 


(93a) 


and 

dq  sinp  sin0  , „ N 

cotp  — - M-de  - f- — r dx  = 0 , (93b) 

^ q cos(0-p)r  ’ 

These  are  the  compatibility  equations  valid  along  the  corresponding 
characteristic  lines.  On  characteristic  lines  the  normal  derivative  of 
the  dependent  variables,  (u,v),  may  be  discontinuous.  This  fact  makes  it 
possible  to  patch  different  flows  together  along  these  lines.  The  only 
restriction  is  that  the  velocity  itself  must  be  continuous.  The  situation 
differs  from  that  for  subsonic  flow,  where  all  derivatives  are  continuous, 
and  a change  in  any  part  of  the  flow  field  affects  every  other  part. 


Along  each  characteristic  line  the  dependent  variables  have  to 
satisfy  the  compatibility  equations.  In  the  two-dimensional  case,  there 
exists  an  analytical  expression  for  the  relation  between  q and  9 but  this 
is  not  true  for  the  axisymmetric  case  and  the  compatibility  equations, 
Equations  (93a)  and  (93b)  have  to  be  solved  in  a step-by-step  procedure 
along  the  characteristic  lines.  The  details  for  calculations  of  the 
interior  and  boundary  points  can  be  found  in  Reference  68  and  are  not 
given  here . 


To  start  the  step-by-step  procedure,  the  flow  properties  p,  T,  M, 
and  0 have  to  be  known  somewhere  in  the  flow  field.  In  the  present  case 
they  are  specified  across  the  exit  plane  of  the  nozzle  along  with  the 
base  pressure,  P^,  the  pressure  distribution  after  the  body,  and  the 

total  temperature,  of  the  exhaust  gas.  Twenty-one  meshpoints  are 

specified  on  the  initial  line.  A concentrated  Prandtl-Meyer  expansion  is 
put  at  the  nozzle  lip  to  match  the  exit  and  base  pressures  and  the  region 
near  the  lip  is  treated  as  a two-dimensional  flow  field.  The  expansion 
fan  is  subdivided  into  10  equal  intervals.  It  has  been  assumed  that  the 
plume  is  axisymmetrical  and  the  calculation  need  only  be  carried  out  for 
the  upper  half  plane.  A consequence  of  this  is  that  configurations  at 
incidence  cannot  be  treated. 

Figure  100  shows  the  effect  of  the  ambient  pressure  changes  on  the 
plume  boundary.  In  one  case,  it  is  assumed  that  the  ambient  pressure  is 
constant  and  equal  to  the  base  pressure,  P^,  while  in  the  other  case  the 
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assumption  was  made  that  the  pressure  recovers  linearly  from  the  base 
pressure  to  the  reattachment  pressure,  P^,  at  the  axial  location, 

x/L  = 0.2.  The  deviations  of  the  resulting  plume  boundaries  are  small 
close  to  the  nozzle  exit  plane. 

Figure  101  indicates  how  the  plume  boundary  expands  as  the  jet  pres- 
sure ratio  is  increased.  When  combined  with  the  strong  interaction  calcu- 
lation procedure,  a new  plume  boundary  is  determined  at  each  step  in  the 
iteration  after  new  values  for  the  base  and  reattachment  pressures  are 
obtained. 

The  present  calculation  is  only  intended  to  provide  the  plume  shape 
close  to  the  body  (at  least  as  far  downstream  as  the  reattachment  point) 
and  no  effort  was  made  to -develop  the  internal  structure  of  the  plume 
in  any  detail.  The  Mach  disc  and  other  shock  situations  arising  inside 
the  jet  are  not  considered.  This  does  not  appear  to  be  a limitation  in 
the  present  application. 

C.  Concerning  the  Nature  of  the  Interaction  Flow 

For  a given  forebody  shape  and  given  nozzle  geometry  it 
is  clear  that  the  flow  in  the  interaction  region  will  be  a function  of 
several  parameters.  Particular  mention  should  be  made  of  the  jet  to 
ambient  pressure  ratio,  the  free  stream  Mach  number,  the  boattail  angle 
and  the  Reynolds  number.  The  influence  of  the  latter  parameter  has  not 
been  greatly  studied  for  the  type  of  applications  of  present  concern. 

It  is  well  documented  that  there  is  a problem  of  great  significance  in 
wing  and  airfoil  design  [48]  and  some  concern  has  been  expressed  in 
relation  to  aircraft  engine  nacelle  and  fuselage  rear  end  design  [67, 

68] . Possibly  for  cases  where  the  jet  pressure  ratio  is  large  there  is 
less  interest  in  a Reynolds  number  effect  and  this  will  not  feature  in 
the  present  discussion;  but  some  theoretical  results  will  be  outlined  in 
a later  section. 

The  experimental  data  obtained  by  Rubin  [ 60]  can  be  displayed  as  on 
Figure  102.  Here,  contours  of  constant  base  pressure  are  contained  on 
a jet  pressure  ratio/boattail  angle  diagram.  The  nozzle  geometry  is 
fixed  (exit  to  base  diameter  ratio  = 0,45;  exit  angle  = 20  ) . The  design 
Mach  number  of  the  nozzle  is  2.7. 

Several  features  of  this  plot  should  be  brought  out.  It  is  seen 
that  each  boattail  configuration  reaches  a minimum  base  pressure  as  the 
thrust  coefficient  of  the  nozzle  is  varied.  In  addition,  the  value  of 
this  minimum  base  pressure  increases  as  the  boattail  angle  is  increased 
but  the  value  of  jet  pressure  at  which  this  condition  is  found  changes 
very  little.  Fortunately,  this  minimum  base  pressure  occurs  at  very  low 
jet  pressure  ratio  and  is  not  of  great  concern  in  the  present  application. 
Physically  the  occurrence  of  this  base  pressure  minimum  implies  that 
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when  the  jet  is  operating  at  low  pressure  ratio  its  influence  on  the 
base  flow  is  similar  to  that  of  an  ejector  and  entrains  more  mass  from 
the  dead-air  region  than  does  the  intrinsic  base  flow  process  (jet-off 
case) . The  additional  mass  entrainment  causes  the  free-shear  layers  to 
turn  through  a larger  angle  and  hence  reduce  the  base  pressure.  The 
effect  of  this  additional  mixing  along  the  outer  edge  of  the  plume  boun- 
dary is  counteracted  as  the  jet  pressure  increases  with  the  consequent 
plume  diameter  increase.  Futher  increase  in  jet  pressure  ratio  results 
in  a monotonic  increase  of  base  pressure  as  the  free-shear  angle  is 
reduced.  But  this  process  cannot  continue  indefinitely.  At  some  point, 
the  free-shear  layer  turns  outward,  instead  of  towards  the  axis  of  the 
flow,  and  the  base  pressure  increases  above  the  ambient  value.  From 
this  point  on,  the  conditions  in  the  free-stream  become  more  significant 
to  the  entire  flow  development. 

Restricting  the  discussion  to  a flow  with  slightly  supersonic  free- 
stream  conditions,  it  is  clear  that  the  outward  turning  of  the  free  shear 
layer  will  sponsor  a shock  wave  in  the  external  flow.  At  this  point  the 
flow  situation  holds  some  analogy  with  that  over  a forward  facing  step 
with  the  step  height  parameter  taking  the  place  of  the  jet  pressure  ratio. 
Czarnecki  and  Jackson  [ 69]  have  presented  experimental  results  for  such 
a flow  and  these  data  demonstrates  that  the  two  flows  do  have  sympathetic 
characteristics . 

For  very  small  step  heights  (corresponding  to  a base  pressure  slightly 
above  ambient),  it  is  found  that  the  external  flow  suffers  an  isentropic 
compression.  For  greater  step  heights,  the  external  flow  compresses 
under  the  classical  shock-boundary-layer  interaction  process.  In  this 
latter  situation  the  pressure  signature  ahead  of  the  step  shows  a rapid 
rise  in  pressure  followed  by  an  essentially  constant  plateau  pressure. 

Returning  to  the  plume  flow,  it  is  expected  that  as  the  jet  pressure 
ratio  increases  to  very  large  values  then  the  base  pressure  would  migrate 
to  a limiting  value  which  associates  with  the  plateau  pressure  for  the 
corresponding  shock-induced  separation  at  the  given  free-stream  Mach 
number  [ 70]  . 

The  base  pressure  contours  in  Figure  102  do  not  reveal  this  limit 
base  pressure  too  clearly  but  the  locus  of  points  at  which  the  external 
flow  first  shows  separation  is  well  defined.  As  may  be  expected,  a lower 
jet  pressure  ratio  is  required  to  separate  the  flow  from  a larger  angled 
boattail  body.  As  was  discussed  in  Reference  19,  this  trend  continues 
for  negative  boattail  angles  (i.e.,  for  flared  bodies).  Figure  104  shows 
the  entire  flow  development  for  a given  body  as  the  jet  pressure  ratio 
is  increased.  Here,  the  free-stream  Mach  number  is  1.1  and  the  boattail 
angle  is  9.5  . The  nozzle  configuration  is  identical  to  that  for  the 
flows  discussed  in  Figure  102.  Both  the  minimum  base  pressure  and  the 
limit  base  pressure  are  displayed  along  with  the  changes  in  flow  structure 
that  ensue  as  the  jet  pressure  ratio  is  increased. 
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The  implication  made  above  that  the  limiting  base  pressure  related 
directly  to  the  plateau  pressure  for  shock-induced  separation  at  the  free- 
stream  Mach  number  needs  some  qualification.  As  the  data  of  Reference 
70  indicated,  this  is  a good  statement  for  a parallel  sided  body.  When 
the  rear  end  is  conical,  however,  the  situation  is  a little  more  complex. 
This  is  due  in  part  to  the  approach  Mach  number  being  different  from  the 
free-stream  Mach  number  and  in  part  to  the  three-dimensional  nature  of 
the  flow. 

Additional  features  of  the  flow  become  evident  from  consideration 
of  Figure  103 . The  free-stream  Mach  number  is  1.1  and  is  close  to  the 
local  Mach  number  ahead  of  the  boattail.  At  low  jet  pressure  ratio,  the 
Mach  number  approaching  the  base  region  will  be  supersonic  and  the  con- 
fluence between  the  free- shear  layer  and  the  jet  is  enacted  under  a shock 
wave.  The  classical  reattachment  analysis  adopted  in  the  above  theoretical 
development  is  essentially  applicable.  As  the  jet  pressure  ratio 
increases,  however,  the  reattachment  (or  wake)  shock  bifurcates  with  the 
forward  limb  moving  onto  the  body  surface  to  play  the  role  of  the  separa- 
tion shock.  The  rearward  limb  passes  on  downstream.  At  these  higher  jet 
pressure  ratio  conditions,  it  is  evident  that  the  confluence  is  no  longer 
taking  place  under  an  oblique  wake  shock  with  supersonic  downstream  con- 
ditions. Most  probably  the  entire  confluence  interaction  is  subsonic  in 
nature.  On  this  basis,  it  is  expected  that  the  theoretical  results  for 
* the  flow  at  high  jet  pressure  ratio  and  low  supersonic  free-stream  Mach 
number  will  be  less  satisfactory. 

Certain  pertinent  pressures  are  a strong  function  of  Mach  number, 
as  has  been  demonstrated  by  the  experimental  work  of  Nash  et  al . [71]  for 
the  transonic  flow  over  a rearward  facing  step.  Figure  104  shows  the 
variation  of  both  the  base  pressure  and  the  reattachment  pressure  with 
Mach  number.  The  peak  pressure  occurring  in  the  flow  is  also  included. 

Very  large  changes  in  the  pressure  levels  of  these  quantities  are  evident 
as  the  Mach  number  passes  through  the  sonic  condition. 

The  three  basic  types  of  interaction  are  shown  diagramatically  on 
Figure  105.  For  high  subsonic  Mach  number  flow  there  is  a slow  increase 
of  pressure  towards  the  separation  point  with  the  entire  interaction 
length  being  quite  large.  The  same  is  true  when  the  approach  Mach  number 
is  at  a low  supersonic  value  and  the  compression  is  isentropic.  These 
two  flows  differ  sharply  from  that  which  occurs  when  the  interaction 
takes  place  under  a shock  wave.  Typical  experimental  data  for  flow  over 
a forward  facing  step  is  shown  in  Figure  106.  Here  it  is  found  that  the 
separation  point  is  very  close  to  the  maximum  pressure  when  the  approach 
conditions  are  subsonic  but  moves  forward  to  adopt  a position  close  to 
the  inflexion  point  in  the  pressure  signature  at  supersonic  conditions. 

The  maximum  pressure  then  takes  over  as  a plateau  value.  Similar  results 
are  found  for  the  surface  pressures  on  a parallel  sided  body  of  revolution 
which  is  suffering  plume  induced  separation  (Figure  107) . These  data  are 
again  taken  from  Reference  60. 
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A correlation  of  the  separation  and  plateau  pressures  was  shown  in 
Figure  93  and  presents  data  for  both  subsonic  and  supersonic  approach 
conditions.  As  seen  from  the  preceeding,  the  subsonic  separation  condi- 
tion passes  over  into  the  plateau  condition  for  supersonic  flow.  The 
supersonic  separation  then  takes  place  at  a relatively  lower  level.  Also 
included  on  that  figure  was  some  data  for  the  limiting  base  pressure  for 
various  boattail configurations . The  influence  of  the  conical  nature  of 
the  flow  is  very  evident  here  for  the  way  in  which  it  dominates  the  level 
of  the  limiting  base  pressure.  All  the  boattail  data  presented  on  Figure 
93  was  for  the  same  nozzle  configuration,  but  it  should  be  recognized 
that  the  nozzle  geometry  does  play  a role  in  the  detailed  development  of 
the  entire  interaction.  A typical  comparison  for  a parallel  sided  body 
is  shown  on  Figure  108.  These  data  shown  in  this  figures  does  not  extend 
to  a sufficiently  high  jet  pressure  ratio  to  determine  the  effect  of 
nozzle  geometry  on  the  limiting  base  pressure. 

Some  of  the  flow  features  discussed  in  this  section  have  considerable 
impact  on  the  development  of  a flow  model;  the  component  model  discussed 
above  in  particular.  The  most  pertinent  of  these  relate  to  conditions 
in  the  separation  and  confluence  regions  of  the  interaction. 

Models  for  flow  separation  under  a small  pressure  gradient  are  not 
well  established.  This  implies  that  the  calculation  procedure  for  the 
entire  interaction  may  not  be  completely  adequate  in  those  situations 
where  the  flow  compresses  isentropieally  to  the  separation  point. 

Typically  this  will  occur  in  that  region  of  the  base  pressure  curve 
(Figure  103)  , where  the  base  pressure  is  just  a little  above  the  ambient 
level.  At  this  stage  in  the  flow  development,  the  compression  waves  are 
coalescing  to  form  the  forward  limb  of  the  wake  shock.  With  further 
increase  of  jet  pressure  ratio,  which  causes  a strengthening  of  the  for- 
ward shock  as  it  moves  ahead  onto  the  body,  this  condition  will  disappear. 

At  the  same  time  as  conditions  near  the  separation  point  become 
less  uncertain  in  their  prediction,  the  confluence  region  deteriorates 
to  a subsonic-like  behavior.  This  difficulty  is  more  -acute  for  bodies 
with  a large  boattail  angle.  In  principle,  the  subsonic  reattachment 
process  is  no  more  difficult  to  discuss  than  is  the  supersonic  one.  The 
problem  is  more  a practical  one  since  appeal  to  oblique  shock  theory  is 
then  no  longer  valid.  It  is  possible  that  the  use  of  the  'small  turning' 
approximation  adopted  herein  is  adequate  for  the  flows  under  discussion 
provided  the  wake  geometry  is  such  that  the  small  disturbance  approxima- 
tions make  sense.  Iteration  between  a numerical  solution  of  the  outer 
potential  flow  and  an  inner  viscous  flow  may  go  some  way  towards  improving 
these  flow  situations  but  a solution  of  the  full  Navier-Stokes  equations 
is  the  long-term  solution.  What  to  do  with  the  Reynolds  stress  tensor 
is  then  the  question  for  discussion. 
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IV.  CALCULATED  RESULTS 


Results  from  the  calculation  procedure  for  plume  induced 
separation  as  outlined  above  have  been  published  in  several  places  [13, 

14,  19,  20,  21,  22,  24]  and  the  significant  features  of  the  calculated 
flow  development  will  be  outlined  below. 

A.  Results  for  Bodies  in  the  Absence  of  Fins 

The  predicted  variation  of  base  pressure  with  pertinent 
parameters  (jet  pressure  ratio,  nozzle  geometry,  etc.)  is  a ready  means 
of  assessing  the  validity  of  the  calculation  procedures. 

As  was  reported  in  Reference  21,  the  nozzle  geometry  is  an  important 
feature  in  the  configuration  design.  These  data  shown  in  Figure  109 
indicate  the  extent  to  which  the  nozzle  geometry  influences  the  varia- 
tion of  base  pressure  with  jet  pressure  ratio  (or  equivalently,  the  thrust 
coefficient) . In  this  comparison  the  nozzle  exit  diameter  ratio  and  the 
nozzle  exit  angle  are  both  varied.  The  nozzle  exit  Mach  number  was 
within  the  range  2.7  - 3.0  and  should  not  constitute  a major  variable  in 
the  comparisons. 

The  effect  of  a single  variable,  the  nozzle  exit  diameter  ratio, 
can  be  isolated  as  shorn  by  the  calculated  results  in  Figure  110.  For 
a jet  pressure  ratio  of  20,  it  is  seen  that  the  base  pressure  is  less 
than  ambient  for  diameter  ratios  less  than  0.9.  The  most  interesting 
feature  of  the  curve  is  the  existence  of  a minimum  base  pressure  when 
the  nozzle  exit  diameter  ratio  is  close  to  0.2.  This  would  appear  to  be 
another  manifestation  of  the  ejector  effect  that  was  discussed  in  con- 
nection with  the  base  pressure  minimum  in  Figure  102.  The  agreement 
shown  between  theory  and  experiment  in  this  situation  is  satisfactory. 

In  both  Figures  109  and  110  the  data  shown  are  for  a parallel  sided 
afterbody. 

The  large  changes  in  base  pressure  that  result  from  changes  in 
nozzle  exit  diameter  ratio  are  reflected  in  the  drag  of  the  body.  The 
theoretical  drag  prediction  is  shown  in  Figure  111  where  the  minimum 
base  pressure  is  returned  as  a maximum  in  the  drag  coefficent  variation. 
The  total  drag  is  also  a strong  and  complex  function  of  the  Reynolds 
number  and  this  subject  will  be  taken  up  later.  In  cases  where  the  rear 
end  is  conical,  the  drag  changes  are  also  more  complex  due  to  the  surface 
pressure  changes  that  are  reflected  as  a boattail  drag.  Again,  details 
of  the  drag  breakdown  will  be  discussed  later. 

The  variation  of  base  pressure  with  increased  jet  pressure  ratio  is 
well  predicted  for  a parallel  sided  body,  as  is  shown  in  Figure  112, 
except  in  the  intermediate  thrust  range  where  the  base  pressure  is  just 
a little  greater  than  the  free-stream  value.  This  relates  to  the  problems 
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discussed  earlier  of  predicting  the  flow  separation  under  a mild  pressure 
gradient  and  the  difficulties  that  present  themselves  when  the  reattach- 
ment process  reduces  to  that  of  a subsonic  flow. 

Figure  113  shows  a comparison  between  the  predicted  and  measured 
base  pressures  for  a body  with  a 6.5  degree  boattail.  This  result  is 
somewhat  disappointing  at  the  higher  jet  pressure  ratios  since  the  theory 
shows  a far  greater  rate  of  attainment  of  the  limiting  base  pressure  than 
is  found  experimentally.  Further  work  is  needed  in  this  area. 

B.  Concerning  the  Forces  on  the  Body 

1.  The  Body  Drag.  The  drag  of  a body  with  a given  geo- 
metric configuration  and  free- stream  Mach  is  a function  of  two  parameters; 
the  Reynolds  number  and  -the  jet  pressure  ratio.  The  Reynolds  number 
affects  the  flow  in  two  ways.  First,  there  i's  the  expected  reduction  in 
skin  friction  coefficient  as  the  Reynolds  number  increases,  and  second, 
there  is  a complex  interaction  with  the  base  flow  due  to  the  changed 
initial  conditions  at  the  start  of  the  interaction.  In  this  connection 
there  is  a strong  cross-influence  with  the  second  parameter-the  jet 
pressure  ratio. 

The  total  drag  is  composed  of  several  components  each  of  which  is  a 
function  of  both  Reynolds  number  and  jet  pressure  ratio.  In  the  present 
theoretical  development,  where  there  is  no  second-order  iteration  between 
the  inviscid  flow  and  the  strong  interaction,  the  form  drag  of  the  fore- 
body is  not  influenced  by  changes  in  either  of  these  parameters  at  0 
incidence.  The  other  form  drag  components  are  influenced  by  variations 
in  both  parameters.  The  skin  friction  is  also  a function  of  the  jet 
pressure  ratio  when  this  is  sufficiently  large  to  move  the  separation 
point  onto  the  body. 

Typical  drag  characteristics  are  presented  on  Figure  114  which 
shows  the  variation  of  total  drag  with  Reynolds  number  for  several 
values  of  jet  pressure  ratio.  For  low  values  of  this  parameter,  it  is 
seen  that  the  total  drag  behavior  is  very  reminiscent  of  that  resulting 
from  the  variation  of  skin  friction  with  Reynolds  number.  At  high  thrust 
coefficient,  however,  the  situation  is  significantly  different  with  the 
total  drag  coefficient  exhibiting  a maximum  at  a Reynolds  number  of  about 
one  million.  As  will  be  seen,  this  drag  reduction  at  low  Reynolds  number 
is  a consequence  of  plume  induced  separation  moving  onto  the  body,  as 
the  initial  conditions  due  to  the  different  boundary  layer  development 
with  Reynolds  number  takes  place.  In  all  of  these  calculations  the  loca- 
tion of  the  transition  point  has  been  held  constant  at  the  shoulder  of 
the  ogive  forebody.  In  a practical  situation,  it  is  expected  that  the 
transition  point  would  also  be  dependent  upon  the  Reynolds  number,  but 
probably  not  too  greatly  due  to  the  favorable  pressure  gradient  along 
the  nose  of  the  body. 

In  a practical  situation  it  is  not  likely  that  the  transition  point 
would  change  its  location  significantly.  It  should  be  noted,  however, 
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that  the  base  pressure  can  be  strongly  influenced  by  changes  in  boundary 
layer  conditions  at  the  start  of  the  interaction.  The  data  presented  in 
Figure  115  shows  the  calculated  base  pressure  for  an  axisymmetric  back 
step  as  a function  of  the  ratio  of  boundary  layer  momentum  thickness  to 
step  height.  The  free- stream  Mach  number  is  again  1.1. 

The  total  drag  distributions  shown  in  Figure  114  are  broken  down 
into  their  several  components  in  Figures  116  and  117.  Figure  117  for  a 
jet  pressure  ratio  of  only  20  is  not  of  great  interest  and  merely  shows 
that  the  total  drag  change  is  mainly  a consequence  of  the  skin  friction 
variation  with  Reynolds  number  under  these  conditions. 

Figure  117  shows  the  breakdown  of  the  total  drag  between  its  various 
components  when  the  jet  pressure  ratio  is  100.  The  situation  is  now  more 
complex  and  the  discussion  will  present  the  flow  development  as  the 
Reynolds  number  is  decreased.  At  this  value  of  jet  pressure  ratio  the 
base  pressure  is  already  close  to  the  separation  value  (being  greater 
than  the  ambient  value,  it  generates  a thrust  on  the  base)  and  increases 
slowly  as  the  Reynolds  number  decreases.  Then,  as  the  Reynolds  number 
reduces  below  1 million,  there  is  a rapid  increase  in  base  pressure  to 
the  plateau  value  as  the  separation  point  moves  forward  onto  the  body. 

The  extent  of  the  separation,  when  present  at  these  lower  Reynolds  numbers, 
is  shown  in  Figure  118.  The  presence  of  the  separated  flow  region  also 
has  an  influece  on  the  skin  friction  drag  (Figure  117) . For  values  of 

6 

Reynolds  number  above  10  there  is  the  expected  increase  of  friction 
coefficient  as  the  Reynolds  number  decreases.  Once  separation  occurs  on 
the  body,  however,  the  forward  movement  of  the  separation  point  is  suffi- 
ciently rapid  to  cause  a reduction  of  the  friction  drag  component. 

The  combined  effect  of  the  forward  movement  of  the  separation  point 
and  the  increased  base  pressure  causes  an  increase  in  the  pressure  level 
over  the  boattail.  This  has  the  effect  of  reducing  the  boattail  drag. 

If  the  Reynolds  number  were  reduced  further  so  that  the  whole  boattail  was 
at  the  plateau  pressure,  then  there  would  be  no  further  ^decrease  of  boat- 
tail  drag. 

These  data  shown  in  Figure  118  is  important  in  its  practical  implica- 
tions: once  the  separation  point  leaves  the  base  of  the  body,  it  moves 
forward  very  rapidly.  This  fact  will  have  severe  consequences  in  relation 
to  the  stability  and  controllability  of  the  body.  An  additional  point 
should  also  be  made  concerning  the  Reynolds  number  variation.  Provided 
the  flow  is  not  too  close  to  the  separation  critical  conditions,  then  it 
is  not  too  sensitive  to  Reynolds  number  (other  than  the  readily  predictable 
change  due  to  skin  friction)  and  no  scaling  problems  between  test  and 
flight  should  be  evident.  However,  if  the  flow  is  in  the  separation 
critical  region,  then  only  a very  small  Reynolds  number  change  can  be 
responsible  for  a vast  change  in  the  flow  configuration.  Then  it  is 

imperative  that  testing  be  carried  out  at  the  correct  Reynolds  number  for 

a realistic  estimate  of  stability  and  controllability  characteristics. 

It  can  be  noted  that  a factor  of  2 in  Reynolds  number  can  cause  a forward 

movement  in  the  separation  point  of  about  157.  of  the  body  length;  the 

same  order  as  the  chord  length  of  the  control  fins. 
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An  additional  comment  concerning  the  drag  of  a body  is  also  of 
interest.  Figures  119  and  120  plot  the  drag  components  against  boattail 
angle  for  a constant  jet  pressure  ratio.  In  general  terms  it  is  found 
that  the  total  drag  is  reduced  as  the  boattail  angle  is  increased.  This 
is  mainly  due  to  the  reduced  pressure  drag  caused  by  the  boattailing.  At 
very  small  boattail  angles,  it  is  the  influence  of  the  boattailing  on  the 
base  drag  that  causes  the  total  drag  reduction  while  at  higher  angles  it 
is  the  boattail  pressure  itself  that  is  responsible  for  the  drag  reduction. 
The  base  drag  changes  for  two  reasons , one  a consequence  of  the  changing 
base  pressure  with  boattail  angle  and  the  other  due  to  the  reduced  base 
area. 


2-.  The  Normal  Force  and  Pitching  Moment.  Figures  121 
and  122  show  the  normal  force  and  pitching  moment  calculated  for  the  MCDAC 
body  (see  Appendix  B)  at  a Mach  number  of  1.0.  These  calculations  are 
for  a body  alone  without' either  fins  or  propulsive  jet. 

As  can  be  expected,  the  agreement  betwen  theory  and  experiment  is 
quite  good  for  both  the  lift  and  pitching  moment  under  these  conditions. 

A more  revealing  test  of  the  method  is  when  a large  jet  pressure  ratio 
plume  is  present. 

No  comparisons  have  been  made  for  the  drag  of  the  configuration  in 
the  present  instance  since  the  test  model  was  sting  mounted  and  the 
computer  program  is  not  complete  for  the  discussion  of  such  flow  situa- 
tions. This  work  will  be  completed  in  the  near  future. 


C.  Other  Considerations 

Besides  such  quantities  as  the  base  pressure  and  forces 
on  the  body,  there  is  also  interest  in  the  length  of  the  interaction  and 
the  size  of  the  plume.  A detailed  study  of  these  has  not  been  made, 
mainly  due  to  the  lack  of  experimental  data.  The  only  available  data 
are  in  the  form  of  Schlieren  photographs  of  which  Figure  123  is  an  example. 
It  is  also  clear  that  the  interaction  geometry  is  not  too  well  defined 
from  such  evidence.  For  comparison,  however,  the  theoretical  estimate 
of  the  flow  shown  in  the  Schlieren  photograph  is  also  included  on  this 
figure.  The  disposition  of  the  wake  shock  does  not  feature  in  the 
theoretical  solution  and  a comparison  of  this  item  is  not  possible. 

The  geometry  of  the  floxtf  in  the  confluence  region  is  shown  in  Figure 
124  for  the  6.3°  boattail  body  and  2 values  of  jet  pressure  ratio.  The 
free- stream  Mach  number  is  1.1.  When  the  jet  pressure  ratio  is  only  50 
then  no  separation  occurs  on  the  body  and  the  confluence  takes  place 
within  one  base  radius  distance  downstream  of  the  nozzle  exit  plane. 

For  the  higher  jet  pressure  ratio  where  separation  is  present,  then  the 
confluence  point  has  passed  further  downstream.  The  length  scale  of  the 
interaction  in  the  streamwise  sense  is  considerably  greater  than  that 
across  the  flow.  This  is  both  in  conformity  with  the  usual  concept  of 
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the  small  turning  approximation  that  was  adopted  in  the  confluence 
region. 

A locus  showing  the  movement  of  the  confluence  point  as  the  jet 
pressure  ratio  is  increased  is  shown  in  Figure  125.  In  this  case  the 
body  has  a 3.2°  boattail.  Of  interest  here  is  the  fact  that  for  small 
jet  pressure  ratios,  there  is  very  little  downstream  movement  of  the 
confluence  point.  The  jet  expands  almost  to  the  base  diameter  before  the 
confluence  point  shows  any  marked  movement  in  the  downstream  direction. 

From  the  interaction  geometries  shown  on  these  figures,  it  is  clear 
that  under  most  situations  the  two  shear  layers  are  in  quite  close  prox- 
imity along  most  of  the  ‘ dead  air"  region.  This  somewhat  reduces  the 
validity  of  the  asymptotic  free  shear  layer  concept  used  in  the  estimation 
of  the  shear  layer  development.  Another  feature  has  been  absent  in  the 
discussion  of  the  base  flow  and  this  concerns  the  effect  of  roll  on  the 
flow.  The  work  of  Reference  74  suggests  that  rotation  of  the  body  can 
have  an  effect  on  the  base  pressure,  due  presumably,  to  the  different 
mixing  rates  that  occur  when  cross  flow  is  present  in  the  shear  layer. 

The  data  from  Reference  74  are  shown  in  Figure  126  and  indicates  the 
effect  of  rotation.  At  low  rotation  speed,  it  is  seen  that  the  base 
pressure  increases,  it  then  reaches  a maximum  and  finally  decreases  rapidly  with 
further  increases  in  rotation  speed. 


Similar  changes  in  the  mixing  rate  of  the  free  shear  layers  could  be 
expected  as  a result  of  the  wake  from  the  fins,  particularly  if  they  have 
a flow  separation.  Since  the  boundary  layer  on  the  fin  has  not  been  dis- 
cussed in  the  present  work,  no  statement  can  be  made  about  the  likely 
effect  of  this. 

D.  Plume  Induced  Separation  of  a Finned  Body 

A few  calculations  have  been  carried  out, to  show  the  effect 
of  fins  on  the  flow  in  the  base  region  and  they  were  reported  in  Reference 
24.  Again,  the  treatment  is  not  rigorous  since  the  flow  was  determined 
along  various  azimuthal  lines  with  the  neglect  of  a cross  flow  component. 

The  procedure  is  as  follows:  the  surface  pressure  distribution  is 
calculated  on  the  body  with  proper  allowance  for  the  effect  of  the  fins 
(within  the  restrictions  of  the  fin  body  theory  used)  and  this  pressure 
distribution  is  used  to  determine  the  plume  interaction  at  a given  azimu- 
thal location.  Thus,  the  presence  of  the  fin  is  included  only  as  far  as 
its  local  effect  on  the  streamwise  development  of  the  boundary  layer,  and 
any  cross  flow  perturbation  is  neglected.  For  small  fin  loading  situa- 
tions, this  may  well  be  an  adequate  treatment  of  the  problem.  It  should 
also  be  noted  that  the  downstream  effect  of  the  fins  is  neglected;  the 
free  shear  layer  mixing  does  not  acknowledge  the  presence  of  the  fins  in 
any  way.  Again  for  small  fin  loading  this  may  not  be  a serious  limitation. 
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Within  these  limitations,  the  effect  of  a fin  configuration  on  base 
pressure  can  be  assessed  from  Figures  127  128.  First,  in  Figure  127, 

is  the  potential  flow  pressures  on  the  base  plane  for  a 2 angle-of-attack, 

with  and  without  fins.  These  base  pressures  are  all  very  close  to  the 
ambient  pressure.  The  effect  of  the  plume  interaction  is  shown  on  Figure 
128  and  it  is  evident  that  the  contribution  from  the  fins  is  minimal 
under  the  conditions  of  this  example. 


V.  CONCLUSIONS  AND  RECOMMENDATIONS 

This  report  has  outlined  the  theoretical  development  of  a 
practical  technique  for  the  estimation  of  the  behavior  of  bodies  of 
revolution  in  transonic  flight.  The  effects  of  both  fins  and  a propul- 
sive jet  have  been  considered  and  results  presented.  It  is  felt  that  as 
an  engineering  tool  the  techniques  laid  down  are  an  adequate  statement  of 
the  state-of-the  art.  Care  should  be  taken,  however,  not  to  construe 
this  remark  as  one  of  complacency.  There  is  still  much  to  be  done  in  the 
area  of  strong  interaction  aerodynamics. 

Thus,  in  several  key  areas  relating  to  the  separation  and  reattach- 
ment process,  appeal  has  been  made  to  empiricism.  This  empiricism  has 
been  based,  in  many  instances,  on  a very  narrow  range  of  data  and  could 
very  well  be  lacking  in  regard  to  general  applicability.  A completely 
numerical  solution  to  the  full  Navier-Stokes  equations  is  the  only  long- 
term endeavor  that  is  likely  to  provide  satisfaction. 
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Figure  1.  Outline  of  flow  configurations 


Figure  2.  Summary  of  appendages  to  the  cylindrical  body  of  revolution 
that  are  addressed  in  the  inviscid  flow  analysis. 
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Figure  5.  Shifting  of  parabolic  point  and  shock 
wave  location  for  various  Mach  numbers. 
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LOCAL  MACH  NUMBER  (Mg) 


Figure  6.  Local  Mach  number  distribution  for  various 

bodies  at  M = 1.0. 
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PRESSURE  COEFFICIENT  (Cp) 


Figure  7.  Effect  of  free-stream  Mach  number  on  the  surface 
pressure  distribution  for  fifth  order  polynomial  body 
and  for  the  ogive-body. 
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PRESSURE  COEFFICIENT  (CJ 
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Figure  8.  Comparison  of  present  transonic  cone  solution 

with  data  6 = 5°,  M = 0.95. 
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PRESSURE  COEFFICIENT  (CJ 


Figure  9.  Comparison  of  present  method 
with  data  for  a blunted  cone. 
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PRESSURE  COEFFICIENT  (Cj 
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Figure  11.  Surface  pressure  distribution  for  a blunted 
cone  at  a 6°  angle-of-attack. 
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Figure  12.  Flow  over  a conical  forebody  6=5° 
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Figure  15.  Comparison  of  theory  with  data  for  flared  body 
with  conical  angle  6 = 6.54°  at  = 0.9. 
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of  theory  with  data  for  flared  body 
at  a = 2°  and  M =0.9. 
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Figure  17.  Variation  of  pressure  coefficient  for  various 

flare  angles  6 at  M =1.0. 
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SURFACE  PRESSURE  COEFFICIENT,  Cp 


Figure  18.  Comparison  of  transonic  similarity  solutions  with  data 
for  15°  half-angle  conical  nose  body. 
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Figure  19.  Comparison  of  zero-order  transonic  similarity  solutions  with  data 
for  10°  half-angle  conical  forebody  at  Mach  number  0.90. 
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Figure  20.  Transonic  similarity  solutions  for 
boattail  with  various  Mach  numbers. 
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Figure  21.  Modified  conical  boattail  program  to  a smooth-shaped 
boattail  for  surface  pressure  distribution. 
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Figure  26.  Computational  model  used  in  Section  II.  D. 
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Figure  28.  Conformal  transformation  for  cruciform  fin-body  combination  [38]. 
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Figure  39.  Comparison  of  load  distribution  on  fin-body 
combinations  in  various  cases. 
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Figure  44.  Total  load  distribution  and  its  elements  (No.  3 fin). 
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Figure  45.  Schematic  model  of  fin-body  combinition  in  present  analysis. 
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Figure  51.  Recalculated  spanwise  lift  distribution 
of  rectangular  fin-body  combination  given  by 
KBrner  [43]. 
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Figure  52.  Longitudinal  pressure  distribution  of  body  with  rectangular  fin. 
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Figure  53.  Chordwise  pressure  distribution  on  45° 
sweptback  fin  with  body. 
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Figure  54.  Spanwise  lift  distribution  with  different 
taper  ratio  fin  given  by  Korner, 
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Figure  56.  Comparison  of  subsonic  similarity  rules  for  cruciform 
slender  delta  fin-body  combination  at  angle  of  incidence  (also 
includes  result  from  slender  body  theory). 
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Figure  57.  Local  lift  coefficient  on  fin  body  junction 
showing  effect  of  Mach  number  change,  -30°  swept  fin 
[43]. 
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Figure  58.  Load  distribution  on  rectangular  wing 
showing  the  effect  of  Mach  number. 
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SLENDER  BODY  THEORY 
(SEE  FIG.  36  OR  REF.  23) 
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Figure  59. 
on  fins 


Comparison  with  slender  body  theory,  distribution  pattern. 
(-6  = 3°,  only  horizontal  fins  are  canted). 


Figure  60.  Comparison  with  slender  body  theory  on  C 
on  body  (-6  = 38,  only  horizontal  fins  are  canted). 
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Figure  61.  Comparison  with  slender  body  theory  of  load  distribution 
pattern  on  fins  (-6  = 3°,  only  horizontal  fins  are  canted). 


(«)  PRESENT  PANELING  SCHEME  (FOR  ZERO  DEFLECTED  FLAP) 


(c)  FLEXSTAB  PANELING  SCHEME.  (Ref.  46) 

Figure  62.  Paneling  schemes. 


© 


CALCULATION) 


2°,  M 

' I 


ffl 


121 


Figure  65.  Flat  wing  with  deflected  flap  (trailing 
edge  6_  =»  8.3°,  a = 0°,  = 0.85). 
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Figure  66.  Normal  force  coefficient 
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Figure  67.  Fin-body  geometry. 
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Figure  68.  Chordwise  C distribution  on  fin:  M “1.62, 
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Figure  70.  Chordwise  load  distribution  on  fin  for  same  conditions 
as  Figure  69.*  Y = 0.752  in. 
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Figure  72.  Chordwise  pressure  distribution  on  fin  for  same 
conditions  as  Figure  69:  Y = 2.555  in. 
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Figure  75.  Comparison  of  compressibility  effect  of  present  result 
with  other  similarity  rules  (a^  = 8.5°). 
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Figure  76.  Variation  of  lift  coefficient  of  fin-body  combination 
with  changing  free-stream  Mach  number. 
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Figure  77.  Diagram  showing  the  major  features  of  the  viscous  flow  interaction. 
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Figure  80.  Calculated  skin  friction  distribution  of  a conical  boattail  for  different 
assumed  Mach  numbers  (M2)  at  the  end  of  the  boattail  shoulder  expansion:  body  B4. 
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Figure  81. 


Diagram  showing  expansion  corner  model  mounted 


in  the  UTSI  transonic  flow  facility. 
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Figure  82.  Measured  surface  pressure  distribution  over  6.5°  expansion 
corner  at  a free-stream  Mach  number  of  0.9. 
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Figure  83.  Measured  velocity  profiles  over  a 6.5°  expansion 
corner  when  the  approach  Mach  number  is  0.7:  corner  is 

located  at  40  in. 
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Figure  84.  Flare  junction  pressure  correlation 
from  the  data  of  Reference  60. 
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Figure  85.  Distribution  of  streamlines  in  shear  layer 
for  various  external  Mach  numbers. 
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DOWNSTREAM 
FLOW  DIRECTION 


Figure  87.  Definition  of  certain  angles  in  the  interaction. 
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Figure  88.  Typical  curve  showing  the  variation  of  base  pressure  with 
the  coefficient  T,  and  the  selection  of  the  empirical  value:  body 

B1  with  nozzle  Nl,  jet  pressure  ratio  = 24.5. 
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Figure  89.  Correlation  of  pressure  rise  to  separation. 
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Figure  90.  Specification  of  the  interaction 
lengths  for  transonic  flow. 
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Figure  91.  Estimates  of  the  plateau  pressure 
for  transonic  flow. 
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Figure  92.  A typical  variation  of  base  pressure  with  jet 
pressure  ratio  showing  the  onset  of  separation  and  the 
limiting  value  of  base  pressure:  body  B1  with  nozzle  Nl. 
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Figure  93(a).  Separation  and  plateau  pressure  ratios  for  supersonic 
flow  over  various  bodies  of  revolution  [60]. 
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Figure  93(b).  Correlation  of  separation  pressure 
for  subsonic  Mach  numbers  - UTSI  data. 
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Figure  94.  Correlation  of  the  limiting  base 
pressure  for  bodies  with  conical  rear  ends. 


Figure  95.  Predicted  pressure  rise  on  a parallel  sided 
body  with  M = 1.1:  body  B1  with  nozzle  Nl. 
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BOUNDARY  LAYER 
DEVELOPMENT, 

SECTIONS  III.A.1  Eon  ,«> 
AND  llt.A.3.  Ad  Em - I55) 
AH  Eqn.  (56) 


Afl,  EON.  (55) 


Figure  96.  Itemized  summary  of  the  calculation  components 
when  the  base  pressure  is  lower  than  ambient. 
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Figure  97.  Summary  of  components  used  in  interaction  calculation 
when  the  base  pressure  is  greater  than  ambient. 
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Figure  100.  Effect  of  various  external  pressure  distributions 
on  plume  boundary  shape. 
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Figure  101.  Effect  of  base  pressure  ratio  on  the  plume  shape. 


0 


158 


BOATTAIL  ANGLE 


Figure  102.,  Contours  of  base  pressure  for  various  boattail 
bodies:  M = 1.1,  bodies  B3,  B4,  and  B5  with  nozzle  N12. 
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Figure  103.  Overall  effect  of  jet  pressure  ratio  on  base  pressure 
and  flow  structure:  = 1.1,  body  B3  with  nozzle  N12. 


Figure  104.  Variation  of  base  pressure  and  reattachment  pressure 
for  the  flow  over  a rearward  facing  step  [ 71] . 
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Figure  105.  Flow  over  a forward  facing  step:  (a)  subsonic 

approach  conditions,  (b)  supersonic  isentropic  flow,  and 
(c)  supersonic  flow  with  shock  separation. 

162 


163 


X/h 

Figure  106(b).  Flow  over  a forward  facing  step  at  = 1.61. 
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Figure  107 . Measured  surface  pressure  distribution 
on  parallel  sided  body:  body  B1  with  nozzle  Nl. 
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Figure  108.  Effect  of  nozzle  geometry  on  the  base 
pressure  signature:  = 1.2,  body  Bl. 
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Figure  109.  The  base  pressure  signatures  corresponding 
to  various  nozzle  configurations. 
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Figure  110.  Calculated  results  for  nozzle  diameter 
ratio  effect  on  base  pressure:  body  Bl. 
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Figure  HI.  Drag  coefficient  variation  corresponding 
to  the  data  in  Figure  110. 
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Figure  112.  Calculated  base  pressure  for  a parallel-sided 
body:  Mw  = 1.1,  body  B1  with  nozzle  N12. 
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Figure  113. 
boattails 


Calculated  base  pressure  on  body  with  a 6.5° 
= 1.1,  body  B4  with  nozzle  N12. 
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Figure  114.  Influence  of  Reynolds  number  on  drag  coefficient 
for  different  jet  pressure  ratios:  M ° 1.1,  body  B4 

with  nozzle  N12. 
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Figure  115.  Effect  of  initial  boundary  layer  thickness 
on  base  pressure  for  flow  over  an  axisymmetric  back 
step:  M = 1.1. 
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Figure  116.  Variation  of  drag  components  with  Reynolds 
number:  = 1.1  and  jet  pressure  ratio  20. 
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Figure  117.  Effect  of  Reynolds  number  on  drag 
components:  Mo  = 1.1  and  jet  .pressure 

ratio  100. 
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Figure  118.  Predicted  separation  point  for  various  jet 

pressure  ratios  as  a function  of  Reynolds  number: 

M = 1.1. 
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DRAG  COEFFICIENT 


Figure  119.  Theoretical  drag  predictions  for  boattailed 
bodies  at  = 1.1:  jet  total  pressure  ratio  is 

50.0;  bodies  B3,  B4,  and  B5  with  nozzle  N12. 
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Figure  120.  Theoretical  drag  predictions  for  boattailed 
bodies  at  = 1.1:  jet  total  pressure  ratio  is  100.0, 
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Figure  123. 
M =1.0: 

CO 


Comparison  of  the  interaction  on  body  B1  with  nozzle  N4 
jet  pressure  ration  = 50.1. 
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Figure  124.  Effect  of  jet  pressure  ratio  on  the  geometry  of  the 
interaction:  body  B4  with  nozzle  N12,  ■ 1.1. 
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Figure  126.  Effect  of  rotation  on  base  pressure  data  [74], 
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Figure  128.  Calculated  b 
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LIST  OF  SYMBOLS 
(Section  II) 

Body  radius,  location  in  x-direction  of  a fin  comer  (Figure 
24). 


Aspect  ratio  of  a fin  extended  into  a body. 

Maximum  spanwise  length,  location  in  y-direction  of  a fin 
corner  (Figure  46) . 

Chordwise  length  at  a fin  root. 

Lift  coefficient. 

Normal  force  coefficient. 


Surface  pressure  coefficient. 


Difference  of  C , i.e.,  AC  = (C  )- 

P ’ ’ p v p' lower 


part 


(C  ) or  (C  ).  - - (C  ) . . „ 

p upper  p left  p right 


part  part 

Body  diameter. 


part. 


Derivative  operator  of  order  m. 

Percentage  of  a location  of  a panel  in  x-direction. 
Location  of  a local  panel  in  y-direction. 

Nonlinear  correction  term. 


Complete  elliptical  integrals  of  the  third  kind  (Reference  23). 


Incomplete  elliptical  integrals  of  the  third  kind  (Reference 
23)  . 

Half  width  of  a local  panel. 


Parabolic  constant. 
Body  total  length. 
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LPC 

i 


M 

M 

CO 

n 

P 

Ap/q 

Q 

q 


R 


1 

Re 


r 

AS 

sra 

n 

s 


U 


00 


W 

X 


Local  panel  chard  length. 

Nose  length  in  II. 9>  chordwise  lengthaty  position  in  II. 3. 
Number  of  panels  on  half  body. 

Mach  number  in  free- stream. 

Number  of  panels  on  half  body. 

Aerodynamic  matrix  [ 45]  . 

Load  distribution  = AC 

P 

Aerodynamic  matrix  £45]  . 

2 

Dynamic  pressure  = p v /2,  source  strength  per  unit  body 
surface  area. 

Body  radius  in  II. 1,  radius  of  circle  in  cr-plane  corresponding 
to  a surface  of  a fin-body  combination  in  x-plane  in  II. 2, 
Aerodynamic  matrix  [45]  in  II. 3. 

Radius  of  secant  body. 

Reynolds  number. 

Polar  coordinate  in  x-plane,  distance  of  two  points. 

Panel  area. 

Sterling  number. 

Semi-spanwise  length  of  a fin. 

Free-stream  velocity. 

Perturbation  velocity  component  in  x-direction. 

Free-stream  velocity.  * 

Perturbation  velocity  component  in  y-direction. 

Velocity  vector. 

Perturbation  velocity  component  in  z-direction. 

Physical  plane  (=  y + iz) . 
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X 


Longitudinal  coordinates . 


£>1,2. 
x,y,z 
x' ,y' ,z' 

Z 

o 

a 

ac 

P 

r 

7 


s 
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0 

0 

« 


9 

9 

p* 

po 

a 

T 

v 


Body -fixed  Cartesian  coordinate  systems  (Figure  22) . 

Fin  location' mounted  on  a body  in  z-direction. 

Angle  of  attack. 

Angle  of  incidence. 

/ 2 

Angle  of  yaw,  ory  1 - M 
Circulation. 

Specific  heat  ratio  = C^/C^,  circulation  strength  per 

unit  spanwise  length,  polar  coordinate  in  a-plane 
corresponding  to  a fin  root  in  x-plane  (Figure  28) . 

Semi-apex  angle  of  nose  cone  in  II. 1,  cant  angle  of  a 
fin  in  II. 2,  deflection  angle  of  a flap  or  elevon  in 
II. 3. 

Polar  coordinate  in  a-plane  corresponding  to  0 in  x- 
plane  (Figure  28) . 

Polar  coordinate  in  x-plane. 

- 90°  -0 

Nonlinear  velocity  potential. 

Linearized  velocity  potential  in  II. 1 or  small  perturba- 
tion velocity  potential  in  II. 2. 

Transformed  perturbation  velocity. 

Angle  of  bank. 

Distance  between  two  points.* 

Density  in  a free- stream. 

Mapping  plane  (=t  + iv) . 

Cartesian  coordinate  in  a-plane. 

Cartesian  coordinate  in  a-plane.  . 
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\ Boattail  angle  in  radians  in  II. 1,  slope  of  a local  bound 

vortex,  taper  ratio. 

A Sweptback  angle  at  a fixed  chordwise  position, 

w Apparent  angle  of  attack  of  a fin. 


SUBSCRIPTS 
b Body, 

f Flap. 

F Fin. 

e Elevon. 

ijjjkjlj'OjtJjp,  panel  name, 
n normal  to  a surface . 

t Tangential  to  a surface, 

d Control  point  of  a panel, 

v Bound  vortes  location. 


a 


erf 

exp 


G 


(Sections  lit  and  IV) 

Sound  speed. 

Drag  coefficient. 

Skin  friction  coefficient. 

Pitching  moment  coefficient. 

Normal  force  coefficient. 

Pressure  coefficient. 

Error  function. 

Exponential  function. 

Modified  trubulence  production  integral. 
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G 

Gcc 

h 

H 

H 

K 

l* 

R(x) 

Re 

8 

T 

T 

u,v 

U 


0 

x.r 


x ,y 
P JP 


M 

P 


P 

q 

y 


Value  of  G for  equilibrium  flow. 

Modified  stream  function. 

Step  height. 

Boundary  layer  shape  factor. 

Incompressible  laminar  flow  value  of  H. 

Constant  (Equation  73) . 

Length  of  free-shear  layer. 

Local  body  radius . 

Reynolds  number  based  on  L. 

Entropy. 

Representation  of  turbulence  stress  (Equation  44) . 
Temperature . 

Velocity  components  along  coordinate  directions  x,r 
Velocity  in  stream  direction. 

Friction  velocity. 

Wake  component  velocity. 

Cylindrical  coordinates. 

Coordinates  of  plume  boundary. 

Length  of  body. 

Mach  number. 

Pressure  gradient  parameter. 

Fluid  pressure. 

Velocity  vector. 

Coordinate  normal  to  shear  layer. 
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7 

6 

5* 

5** 

r 

p 

x 

e 

e* 

6M 

P 

♦ 

O 

T 

X 

it 


Non-dimensional  interaction  length  scale. 

Ration  of  specific  heats. 

Boundary  layer  thickness . 

Displacement  thickness  of  boundary  layer. 

Density  thichness. 

Pressure  rise  coefficient  (Equation  69). 

Fluid  density. 

Non- equilibrium  parameter  or  characteristic  angle. 
Boundary  layer  momentum  thichness  or  turning  angle. 
Boundary  layer  energy  thickness. 

Angle  for  shock  induced  separation 
Mach  angle. 

Stream  function. 

Kinematic  viscosity  of  fluid. 

Reynolds  shear  stress. 

Von  Karman  constant. 

3.14... 


A Modified  shape  factor. 

T)  Similarity  parameter. 

a Shear  layer  mixing  cofficient. 

6 

X Pressure  signature  function. 

SUBSCRIPTS 

B Base  condition. 

BL  Boundary  layer  . 

C Chamber  condition  or  confluence  point  . 
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e 

f 

J 

m 

o 

P 

ref 

S 

SL  or  S 

t 

w 

CO 

2 


Condition  external  to  boundary  layer. 
Condition  at  flare  juction. 

Jet  value. 

Value  on  median  streamline. 
Stagnation  or  initial  condtion. 
Plateau  or  plume  value. 

Reference  condition. 

Separation  condition. 

Free  shear  layer. 

Value  at  transition. 

Value  on  wall. 

Free  Stream  Condition. 

Value  downstream  of  discontinuity. 
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Appendix  A.  COMPUTER  PROGRAM  UTILIZATION 


A subsequent  document  will  give  greater  details  of  the  computer 
program,  along  with  a full  listing  and  documentation  of  its  capabilities. 
The  objective  here  is  to  simply  outline  the  input  to  the  program  and  to 
present  a sample  output.  Table  A-l  gives  the  symbol  and  format  inputs 
for  calculation  of  the  potential  flow  and  the  viscous  interaction.  Table 
A-2  presents  the  same  information  for  the  calculations  when  the  surface 
Mach  number  distribution  is  given.  In  this  latter  case,  there  are  again 
two  options  [with  IOPT  being  redefined  (see  Table (A-2) ] so  that  either  a 
boundary  layer  calculation  or  a combined  boundary  layer  calculation/plume 
interaction  can  be  determined.  This  option  allows  experimental  pressure 
distributions  along  the  body  to  be  used  in  place  of  a potential  flow 
calculation.  The  input  surface  pressure  distribution  can  be  in  the  form 
of  a pressure  coefficient,  when  the  parameter  IM  = 1,  or  as  a local  Mach 
number  distribution  (then  IM  = 0) . All  internal  calculations  in  the 
program  use  the  local  Mach  number  and  all  required  interchanges  between 
Mach  number  and  pressure,  or  velocity  are  done  to  the  basis  of  isentropic 
flow. 


The  basic  flow  diagram  through  the  program  is  shown  in  Figure  A-l, 
while  Table  A-3  explains  the  notation  used  in  Table  A-l  and  A-2  and  in 
the  computer  program.  Reference  should  also  be  made  to  Figure  A-2. 

All  geometric  lengths  associated  with  the  body  configuration  are 
input  in  the  same  units  and  non-dimensionalized  (with  respect  to  BDLN) 
within  the  program.  Non-dimensional  input  is  admissible  provided  BDLN 
is  specified  to  be  unit  length.  Angles  are  input  in  degrees. 

The  potential  flow  calculation  requires  the  specification  of  two 
parameters  CS  and  A.  Appropriate  values  for  these  are  given  in  Table 
A-4.  During  the  process  of  the  calculation,  there  is  an  iteration  on 
these  values  to  determine  the  final  solution.  The  final  values  for  CS 
and  A are  displaced  (Figure  A-3a)  where  CS  identifies  with  the  sonic 
point  and  A retains  its  symbol.  As  experience  is  gained  in  the  use  of 
the  program  for  various  body  shapes,  better  first  estimates  for  CS  and 
A than  those  provided  in  Table  A-4  can  be  established  for  faster 
convergence . 

Sample  data  cards  are  shown  in  Table  A-5  and  the  corresponding 
computer  output  is  given  on  Figure  A-3.  In* this  example  there  are  four 
figures  of  output  giving  the  potential  flow,  the  boundary  layer  develop- 
ment and  the  interaction  with  the  plume  (Figure  A-3b-A-3e) . The  body 
selected  for  this  example  had  a 4- calibre  tangent  ogive  nose  and  a 
parallel  sided  aft-end.  There  were  no  fins.  It  should  be  noted  that 
the  geometric  data  of  the  body  are  returned  to  the  output  in  non- 
dimensional  form.  The  body  length  (BDLN)  is,  however,  recorded  in  the 
output.  Notation  used  with  the  output  display  is  presented  on  Table  A-6. 
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The  option  IC  = 0 or  1 determines  whether  the  potential  flow  is 
calculated  or  not.  The  potential  flow  is  not  calculated  when  a surface 
pressure  distribution  is  input  data  (I0PT  = 3) . The  main  value  of  the 
XC  = 0 option  is  that  it  allows  several  (ICASE)  different  viscous  flow 
calculations  (with  various  Reynolds  number  or  jet  pressure  ratio)  for  only 
one  potential  flow  calculation.  Thus,  put  IG  = 1 for  the  first  data  set 
and  the  specification  IC  = 0 would  be  made  for  all  subsequent  data  sets 
in  the  sequence.  The  potential  flow  solution  is  stored  in  the  program 
until  overwritten  by  a subsequent  calculation  or  input  data  in  the  I0PT  = 

3 mode. 

The  program  is  normally  run  with  the  calculation  being  performed 
along  a single  azimuthal  strip  at  the  specified  angle  BE.  In  this  mode 
IL  = 0,  Placing  IL  = 1. in  the  input  data  will  initiate  a sequence  of  cal- 
culations for  BE  = 0,  30,  60,  ...,  180  and  an-  integration  for  the  normal 
force  coefficient. 


TABLE  A-l.  INPUT  CARDS  FOR  POTENTIAL  FLOW  AND  VISCOUS  INTERACTION 
• CALCULATIONS  (I0PT  =0,  1 OR  2) 


Card 

Symbol  Sequence 

Format  (Fortran  IV) 

ICASE 

13 

TITLE  CARD  FOR  EACH  RUN* 

18A4 

I0PT,  IC,  IF,  IND,  NST,  IL 

615 

MINF,  CS,  A,  PINF,  REN 

5F  10.6 

BDLN,  XSHLD , XBT,  XN,  XS 

5F  10.6 

REN0Z,  PJ,  THJ,  XMJ-,  THJl,  T0 

6F  10.6 

THBL,  THC,  RB,  RC 

4F  10.6 

8 

AL,  BE,  THPS , RCN 

4F  10.6 

9 

ANG,  CRD,  FSPN,  XF,  PH 

5F  10.6 

* Cards  2-9  are  repeated  "ICASE"  times 


TABLE  A- 2.  INPUT  CARDS  FOR  BOUNDARY  LAYER  AND  VISCOUS  INTERACTION 
CALCULATIONS  (I0PT  = 3) 


Card 

Symbol  Sequence 

Format  (Fortran  IV) 

1. 

ICASE* 

13 

2-9 

As  in  Table  A-l 

As  in  Table  A-l 

10 

IN,  IM,  I0PT** 

313 

11 

XST 

F 10.6 

12 

X(I)  (1=1,  IN) 

7F  10.6 

13 

R(I)  (1=1,  IN) 

7F  10.6 

14 

M(I)  (1=1,  IN) 

7F  10.6 

*Cards  2-14  are  repeated  "ICASE"  times 

**  Only  values  I0PT  = 0 or  2 are  sensible 
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TABLE  A-3,  NOTATION* 


Symbol 

Definition 

Units 

ICASE 

Number  of  runs  in  given  data  set. 

- 

10  FT 

Selects  calculation  options  (See  Figure 
A-l) . 

~ 

IC 

IC  = 0 omits  the  potential  flow 
calculation. 

- 

IF 

IF  = 0 omits  the  fin  calculation. 

- 

IND 

IND  = 1'  for  potential  flow  calculations 
on  the  body. 

\ 

IND  = 2 is  for  calculations  off  the  body 
(done  at  radius  RCN)  but  within  a mixed 
subsonic/ supersonic  region. 

IND  = 3 does  off  body  calculations  in 
cases  where  the  flow  character  does  not 
change . 

NST 

NST  = 0 for  base  flows  with  a plume. 

- 

NST  = 1 for  base  flow  in  the  absence  of 
a plume. 

- 

IL 

IL  = 0 calculates  flow  at  one  azimuthal 
angle  only. 

- 

IL  = 1 determines  flow  at  seven  azi- 
muthal angles  and  determines  the  lift 
coefficient  (a  * 0 only) . 

MINF 

Free- stream  Mach  number. 

- 

CS,  A 

Constants  for  potential  flow  calculation 
(See  Table  A-4) . 

PINF 

Free-stream  pressure. 

psi 

REN 

Reynolds  number  in  millions. 

- 

BDLN 

Reference  length  of  body. 

Length 

(units  L say) 

* See  also  Figure  A-2. 
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TABLE  A-3 . (CONTINUED) 


Symbol 

Definition 

Units 

XBT 

Distance  from  nose  to  boattail  (or 
flare)  corner. 

L 

XN 

Nozzle  length. 

L 

XS 

Transition  point  for  boundary  layer 
calculation. 

L 

REN0Z 

Nozzle  exit  radius. 

L 

PJ 

Jet  total  pressure  ratio. 

- 

TRJ 

Nozzle  exit  angle. 

Degrees 

XMJ 

Nozzle  exit  Mach  number. 

- 

THJl 

Boundary  layer  momentum  thickness  at 
nozzle  exit. 

L 

T0 

Jet  total  temperature. 

°R 

THBL 

Boattail  conical  angle  (<  0 for  flare) . 

Degrees 

THC 

Nose  conical  angle  (conical  nose  only) . 

Degrees 

RB 

Base  radius. 

L 

RC 

Main  body  radius. 

L 

AL 

Bo(dy  incidence. 

Degrees 

BE 

Azimuthal  angle  (see  Figure  4) . 

Degrees 

THPS 

Initial  shear  layer  angle.  . 

Radians 

RCN 

Off-body  radius. 

L 

ANG 

Fin  cant  angle. 

Degrees 

CRD 

Fin  root  chord. 

L 

FSPN 

Fin  span. 

L 

XF 

Location  of  fin. 

L 

PH 

Azimuthal  location  of  fin  (see  Figure 
22)  . 

Degrees 

IN 

Array  size  for  input  pressure 
distribution. 

IM 

IM  = 0 if  Mach  number  distribution 
input . 

- 

IM  = 1 if  pressure  coefficient  distri- 
bution input. 

- 
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TABLE  A-3 . (CONCLUDED) 


Symbol 

Definition 

Units 

XST 

Location  of  discontinuity  in  input 

pressure,  distribution. 

L 

X(I) 

Distance  array  for  input  pressure 

L 

distribution. 

(I  = 1,  IN) 

R(I) 

Body  radius  array  for  input  pressure 

distribution. 

(I  = 1,  IN) 

M(I) 

Input  pressure  distribution  (see 

definition  for  IM) 

(X  = 1,  IN) 

TABLE  A-4.  VALUES  FOR  CONSTANTS  CS  AND  A 


MINF 

CS 

A 

MINF  < 1.0 

0.67 

-17.5 

MINF  =1.0 

0.5 

- 2.0 

MINF  >1.0 

0.5 

-10.0 
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TABLE  A-6.  SYMBOLS  USED  ON  COMPUTER  OUTPUT 


Symbol 

Definition 

X 

Non-dimensional  distance  along  body. 

X (page  4 only) 

Non-dimensional  distance  along  plume 
(origin  at  nozzle  exit  plane) . 

R 

Non-dimensional  local  body  radius. 

R (page  4 only) 

Non-dimensional  plume  radius. 

M 

Local  Mach  number. 

DR/DX 

Body  slope. 

CP 

Local  pressure  coefficient. 

D 

Body  reference  parameter  (=  2 RC) . 

P 

Local  pressure  (psi) . 

PINE 

Free-stream  pressure. 

H 

Boundary  layer  form  factor. 

THETA 

Boundary  layer  momentum  thickness  (non- 
dimensional)  . 

DEL 

Boundary  layer  displacement  thickness 
(.non-dimensional)  . 

CF 

Local  skin  friction  coefficient. 

BLTH 

Boundary  layer  thickness  (non-dimensional). 

S 

Non-dimensional  distance  along  plume 
boundary . 

XP,  YP 

Coordinates  of  confluence  point  (non- 
dimensional)  with  origin  at  nozzle  exit 
plane. 

209 


210 


Figure  A-l.  Basic  flow  diagram  for  computer  program. 
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Figure  A -3 a. 
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Figure  A-3b.  Boundary  layer  calculation. 
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Figure  A-3c 

. Pressure  distribution  with  plume  interaction. 
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SAMPLE  CALCULATION 

^FREE  STREAM  MACH  hUHkEp  = 1 « x 0 0 0 
JFT  PPFSSURF  RAT lu  = .,.7667 

JFf  ‘ TOTAL  PRESSURE  - 67.7000 

THRUST  COEFFICIENT  = Ot'.Er.rio 
SEPARATION  POINT  AT  l.uOOU 
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free’  stream  PRESSURE  =•  14.7000 
JFT  TOTAL  TEMPERATURE  = 500.00 

BASE  PRESSURE  = 1.277a 

JET  EXIT  MACH  NUMBER  = 2.7000 
PLUME  OtUMETRY 


t 


l 

1 . 

X 

1 

1 

1 

R 

1 

1 

1 

S 

1 

1 

1 

1 

1 

1 

1 

1 

0 . 000000 

1 

. o.'isc)  11 

1 

0.000000 

1 

i 

• 0i>39a6 

1 

• 03.1464 

1 

.004744 

1 

■1 

.012295 

1 

. 03-446 

I 

.014492 

1 

1 

. 0 i 645H 

1 

. 04ij573 

1 

.019144 

1 

l 

•02072a 

1 

.048507 

1 

.023828 

1 

1 

• 0 251 X 1 

1 

.044315 

1 

.028573 

I 

I 

.029^,21 

1 

.04-S999 

1 

.033387 

1 

i 

.034254 

J 

. .047559 

1 

.038276 

1 

l 

.039010 

1 

. 0AH992 

l 

.043243 

1 

X 

.<}438m7 

1 

. 051-295 

1 

.048291 

.1 

i 

• U468d4 

1 

.051465 

I 

.053423 

1 

1 

*u53998 

1 

• 052446 

1 

.058640 

I 

i 

.059225 

1 

- .055384 

1 

.063943 

1 

i 

« 0o45o.3 

1 

.054121 

1 

.069330 

1 

X 

.070005  . 

1 

.054700 

1 

.074803 

X 

1 

•075546 

1 

.0551  IS 

1 

.080360 

1 

l 

.OB 11 80 

1 

.055356 

k 

' .085999 

1 

i 

.086899 

1 

.056414 

1 

.091719 

1 

X 

.092696, 

1 

.056281 4 

1 

.097517 

X 

1 

.09B5o2 

1 

.054946 

1 

.103392 

1 

1 

1 

- 

1 

1 

NOZZLE  HALF"jui6Lt  =23.0uo0 
NOZZLE,  RADIUS  = .0306 

BASE ■ RADIUS  = .0J8S  .DIAMETER  RATIO  = ' .8016 

PLUME  INTERSECTION  point  XP  = 1,0158  YP  = »0A02 

REATTACHMENT  PRESSURE  = l.j-,3 

PLATEAU  PRESSURE  = 1.4-.9J 
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Figure  A-3d.  Details  of  the  interaction. 
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sample  calculation 


Summary  of  body  geometry 


BODY  LENGTH  -L  =32.5000 
BODY  R ADIUS/L  = . .0355 
SHOULOER  LOCATION/L  = .3077 

NOZZLE  EXIT  PLANF/L  ='  1.0000 


fHEE  STREAM  CONDITIONS 


INCIDENCE  = 0.0000  DEGREES 
FREE  STREAM  “ACM  NUMBER  = I. 1000 
REYNOLDS  NUMBER  = 1.2000  MILLION 
JET  EXIT  MACH  NUMBER  = 2.7000 
JET  TOTAL  PRESSURE  RATIO  =2A.5000 


SUMMARY. OF  FORCE  COEFFICIENTS 


normal  force  coefficient  = 0.0000 

PITCHING  MOMENT  COEFFICIENT  ABOUT  NOSE  = 0.0000 
TOTAL  DRAG  COEFFICIENT  = .2721 


DRAG  COMPONENTS 

PRESSURE  FRICTION 

BOATtAIL 

BASE 

eno  of  calculation 

.0562  .1272 

O.OaOO 

.0947 

Figure  A-3e.  Summary  of  the  flow. 
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Appendix  B.  SUMMARY  OF  BODY  GEOMETRIES 
USED  IN  STUDY 


All  bodies  "B"  are  fitted  with  a 4-caliber  tangent  ogive  nose  and 
have  a 2.5-in.  diameter  cylindrical  main  body.  The  entire  length  of  the 
bodies  was  32.5  in.  MCDAC  is  a 3-caliber  tangent  ogive  nose  body  which 
is  37.5-in.  long  and  3.75  in.  in  diameter.  A summary  of  bodies  are 
presented  in  Table  B-l,  and  a summary  of  nozzles  are  presented  in  Table 
B-2. 


TABLE  B-l.  BODIES 


Symbol 

Boattail  Angle* 
(deg) 

Boattail  Length 
(in.) 

Bl 

0.0 

B3 

9.5 

2.05 

B4 

6.3 

2.05 

B5 

3.2 

2.05 

B7 

-3.3 

7.50 

B8 

-6.5 

3.75 

B9 

-3.3 

3.75 

MCDAC 

0.0 

- 

*Flared  bodies  are  donated  by  a negative 
angle. 


TABLE  B-2.  NOZZLES 


Symbol 

Diameter  Ratio* 

Exit  Angle, 
(deg) 

Exit  Mach  Number 

Nl 

0.800 

20 

2.7 

N4 

0.78 

0 

2.7 

N12 

0.450  • 

20 

2.7 

N14 

0.600 

0 

2.7 

*Based  on  2,5  in.  cylinder  body  diameter 
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